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APPLICATION  OF  NONDESTRUCTIVE  TESTING  TECHNIQUES  TO  MATERIALS  TESTING 

G.  S.  Kino 


INTRODUCTION 

During  the  last  year,  we  have  demonstrated  a  range  of  fundamentally  new 
optical  nondestructive  techniques  for  measurements  of  surface  defects,  range, 
surface  profiles,  and  acoustic  and  thermal  displacements  of  a  surface.  We 
believe  that  these  will  be  of  great  Importance  both  for  nondestructive  testing 
and  for  manufacturing  technology. 

During  our  initial  research  on  this  contract  on  the  scanning  acoustic  mi¬ 
croscope,  we  developed  several  new  quantitative  measurement  techniques.  He 
have  now  shown  these  same  techniques  can  be  applied  to  optical  microscopy  and 
to  optical  range  sensing.  Mdch  to  our  surprise,  coherent  optics  has  not  been 
used  In  this  way  before.  We  have  therefore  decided  to  exploit  the  advantages 
of  optics.  Its  relative  insensitivity  to  environmental  conditions,  its  very 
good  resolution,  and  no  contact  with  the  sample.  We  also  plan  to  develop  a 
new  range  of  tools  for  range  sensing,  profile  measurements,  flaw  detection, 
and  nondestructive  testing  where  the  optical  probe  Is  used  as  a  thermal  wave 
or  acoustic  wave  sensor. 

In  our  acoustic  microscope  work,  we  have  demonstrated  that,  by  using  a 
wide  aperture  acoustic  beam,  we  can  measure  range  to  a  far  better  accuracy 
than  with  conventional  NDT  systems  using  short  pulses.  The  range  resolution 
of  such  a  wide  aperture  scanning  microscope  lens  Is  comparable  to  a  wave¬ 


length;  the  range  resolution  of  a  short  pulse  system  Is,  at  best,  two  or  three 
wavelengths.  This  fundamental  breakthrough  made  It  possible.  In  the  work  of 


Khurl-Yakub,*  to  measure  the  properties  of  near-surface  flaws  without  In¬ 
terference  from  surface  reflections.  Khurl-Yakub  has  exploited  these  tech¬ 
niques  for  NOT  of  composites  and  of  near  surface  flaws  In  metals,  glass,  and 
ceraml cs . 

Normally,  optical  systems  do  not  give  very  good  range  resolution,  but  the 
use  of  this  same  phenomenon  has  made  It  possible  to  develop  new  optical  micro¬ 
scopy  tools  for  measuring  the  thickness  of  thin  films.  Since  we  can  also  mea¬ 
sure  amplitude  and  phase  of  such  coherent  waves,  we  can  In  fact  measure  the 
distance  of  a  surface  from  the  lens  with  accuracies  of  the  order  of  .001  x  . 

In  our  acoustic  microscope  work,  operating  at  a  wavelength  of  30  ym  ,  we  were 

able  to  measure  film  thicknesses  to  an  accuracy  of  an  order  of  200  A  .  In 
our  optical  microscope  work,  working  at  a  wavelength  of  5100  A  ,  we  are  able 
to  measure  film  thicknesses  to  an  accuracy  of  a  few  Angstroms. 

These  basic  developments  have  led  to  a  seed  program  supported  by  SIMA, 
the  Stanford  Institute  for  Manufacturing  and  Automation,  for  measurements  on  a 
much  larger  scale.  He  have  demonstrated  a  range  sensor  for  robotics  with  a 
range  accuracy  of  the  order  of  0.5  mm  at  20  cm  ,  and  another  type  of  more 
accurate  range  sensor  with  accuracies  of  the  order  of  a  few  microns  at  dist¬ 
ances  of  10-20  cm  .  This  latter  type  of  range  sensor  should  be  of  great 

Importance  In  applications  to  precision  machining.  We  now  envisage  the  possi¬ 

bility  of  carrying  out  measurements  with  the  accuracy  of  a  high-quality  mi¬ 
crometer  by  using  noncontacting,  nondestructive  optical  methods.  Furthermore, 
the  fundamental  tools  that  we  are  developing  will  not  only  enable  us  to  mea¬ 
sure  range  directly,  but  also  to  measure  such  parameters  as  the  slope  of  a 

*P.  Relnholdtsen  and  B.  T.  Khurl-Yakub,  "The  Effects  of  Surface  Roughness  on 
Subsurface  Defect  Detection  Using  Acoustic  Microscopy,"  Presented  at  the 
DARPA/AF  Review  of  Progress  In  Quantitative  NDE,  Williamsburg,  VA  (July  1985). 


surface,  surface  roughness,  and  surface  curvature,  all  without  making  contact 
with  the  sample. 

A  few  examples  may  point  out  the  importance  of  these  possible  break¬ 
throughs  In  nondestructive  testing  technology.  Consider  the  problem  of  mea¬ 
suring  the  Inside  of  an  engine  cylinder  during  manufacture  and  at  a  later 
stage  when  testing  It  In  the  field.  We  envisage  making  optical  measurements 
of  the  cylinder  radius,  how  much  out  of  round  It  Is,  Its  surface  finish,  and 
whether  there  are  flaws  In  the  surface.  We  also  envisage  carrying  out  mea¬ 
surements  during  the  machining  process  Itself  without  touching  the  object. 

The  same  technique  could  also  be  used  to  make  profile  measurements  of  large 
structures  like  an  automobile  or  a  complicated  engine  component.  On  the  other 
hand,  this  kind  of  profiling  measurement  can  be  made  on  a  small  scale  on  semi¬ 
conductors  and  Integrated  circuits. 

Most  of  our  concentration  In  this  work  has  been  on  techniques  suitable 
for  measuring  small  electronic  components,  precision  mirrors,  or  magnetic  re¬ 
cording  disks.  The  basic  principles  we  are  developing,  however,  apply  both  to 
small  and  large  structures,  and  we  have  carried  out  enough  demonstrations  to 
show  that  the  principles  are  universal. 

We  use  laser  sources  for  the  purpose  because  of  their  high  Intensity.  In 
many  of  the  techniques  we  have  developed,  however,  incoherent  light  of  high 
enough  Intensity  would  do  as  well.  The  Implication  of  this  statement  Is  that 
the  surfaces  we  are  looking  at  need  not  be  mirror  smooth.  In  the  larger  scale 
applications,  the  techniques  should  work  with  machined  surfaces  and  materials 
which  are  not  necessarily  lapped  to  a  fine  finish.  We  have  made  basic  demon¬ 
strations  to  show  that  these  conclusions  are  valid. 

Thus,  In  this  work,  we  are  developing  a  range  of  new  techniques  based  on 
our  original  work  with  the  acoustic  microscope,  which  have  now  proven  to  be 
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very  useful  In  optical  applications.  In  these  applications.  Me  can  measure 
range,  curvature,  electronic  properties,  and  material  properties,  and  we  can 
scan  an  optical  beam  at  will.  The  mere  fact  that  we  can  scan  the  optical  beam 
and  produce  many  optical  spots  at  once  also  leads  to  possible  new  techniques 
for  photolithography  and  for  semiconductor  processing  with  optical  beams.  The 
range  of  possible  applications  Is  very  large,  but  they  are  based  on  the  devel¬ 
opment,  In  every  case,  of  nondestructive  techniques  which  do  not  need  any 
contact  with  the  sample  and  which  can  be  used  at  high  speed. 

We  believe  that  these  techniques  are  major  breakthroughs  In  the  field  of 
nondestructive  testing. 

ACOUSTIC  MICROSCOPY 

This  work  Is  now  finished.  Three  papers  published  on  the  work  are  given 
as  Appendix  I  (332,333  C-H  paper)  In  this  report.  They  describe  the  results 
that  we  obtained. 

There  were  two  basic  developments  In  this  research.  The  first  was,  for 
the  first  time,  to  measure  both  the  amplitude  and  phase  of  the  received  signal 
In  an  acoustic  microscope.  Since  we  could  measure  phase  to  a  tenth  of  a  de¬ 
gree,  this  made  It  possible  to  measure  heights  to  an  accuracy  of  a  few  thou¬ 
sandths  of  a  wavelength.  Thus,  we  were  able  to  measure  film  thicknesses  to 
accuracies  of  the  order  of  200  A  using  an  acoustic  wavelength  of  30  pm  In 
water.  More  Important  still,  for  the  first  time  we  were  able  to  make  stress 
measurements  In  glass  samples  with  a  resolution  of  the  order  of  one  or  two 
millimeters  over  the  surface  of  a  tempered  glass  sample.  This  technique  may 
prove  to  be  of  great  Importance  for  measuring  near-surface  stress  variations 
In  ceramics  and  other  materials  In  which  there  are  Inhomogeneities  present. 


The  second  major  development  was  based  on  the  use  of  the  so-called  V(z) 
measurement.  As  the  lens  of  an  acoustic  microscope  Is  moved  towards  the  sam¬ 
ple  being  examined,  the  focus  moves  from  the  point  above  the  surface  of  the 
substrate  to  a  virtual  point  below  the  surface.  If  the  distance  from  the  sur¬ 
face  to  the  focus  Is  called  z  ,  the  return  signal  amplitude  at  the  detector 
Is  called  V(z)  .  If  we  plot  V(z)  as  a  function  of  z  ,  we  typically  obtain 
a  quasi -periodic  curve  with  the  distance  between  minima  dependent  on  the  Ray¬ 
leigh  wave  velocity  of  waves  along  the  surface  of  the  substrate.  We  developed 
a  theory  which  showed:  (1)  the  relation  between  V(z)  and  the  reflectivity  of 
the  sample  as  a  function  of  angle  of  Incidence  of  a  plane  wave,  and  (2)  that 
the  V ( z )  plot  could  be  Inverted  to  determine  the  reflection  function  of  a 
plane  wave  from  a  surface. 

These  results  made  It  possible  to  use  an  acoustic  microscope  to  directly 
measure  not  only  the  surface  wave  velocity,  but  the  shear  wave  and  longitud¬ 
inal  wave  velocities  of  waves  along  the  surface.  Thus,  we  could  use  the  mi¬ 
croscope  to  Investigate  the  material  properties  In  a  small  region  and  put  the 
results  on  a  quantitative  basis.  The  theory  and  experiments  were  In  excellent 
agreement. 

The  results  were  very  well  received  by  the  acoustic  microscope  community, 
for  this  was  the  first  demonstration  that  a  complete  Inversion  of  the  V(z) 
plots  could  be  obtained.  Now  the  theory  Is  being  further  developed  for  use 
with  optical  microscopy  and  range  sensing. 

OPTICAL  MEASURE  f€NTS 
Introduction 

The  scanning  acoustic  microscope  used  coherent  acoustic  waves  to  produce 
Images  free  of  speckle,  and  Is  capable  of  quantitative  phase  and  amplitude 
measurements.  These  results  have  led  us  to  Investigate  how  these  techniques 
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might  be  employed  In  optics.  Interestingly  enough,  quantitative  optical  mea¬ 
surement  techniques  of  this  type  have  not  been  carried  out  before.  The  optic¬ 
al  projects  on  this  contract  have  therefore  been  aimed  at  investigating  how  we 
could  Improve  optical  systems  by  employing  modern  electronic  techniques  within 
the  optical  system.  Our  purpose  has  been  to  make  accurate  profile  measure¬ 
ments  with  better  transverse  resolution  than  Is  presently  obtainable  with  the 
optical  microscope,  and  to  make  quantitative  measurements  of  height  rather 
than  just  produce  a  qualitative  phase  contrast  Image. 

We  demonstrated  an  electronically-scanned  optical  microscope  which  could, 
for  the  first  time,  measure  both  optical  phase  and  amplitude  simultaneously. 

We  used  these  results  to  obtain  vertical  resolutions  of  the  order  of  10  A  or 
better  and  profiles  with  a  transverse  resolution  of  0.27  ym  .  Because  we  had 
optical  phase  and  amplitude  Information  available,  we  were  able  to  demonstrate 
that,  by  Fourier  transform  processing,  we  could  Improve  the  transverse  resolu¬ 
tion  to  0.13  um  . 

These  results  were  an  Important  breakthrough  because  Improvement  In  the 
resolution  of  the  optical  microscope  and  Its  accuracy  for  quantitative  mea¬ 
surements  are  badly  needed  for  measuring  small  components,  such  as  those  used 
In  electronics.  They  are  also  needed  for  precision  machining  applications, 
such  as  the  machining  of  optical  mirrors  and  the  determination  of  small  sur¬ 
face  defects  In  a  wide  range  of  precision  machined  parts. 

With  these  early  results  In  hand,  we  were  then  able  to  go  on  and  incor¬ 
porate  other  electronic  components  In  an  optical  system  which  could  be  used 
for  optical  range  sensing,  for  slope  measurements ,  and  for  edge  measurements. 
It  would  also  be  feasible  to  use  the  system  for  measurements  of  curvature  and 
other  parameters.  Later,  we  realized  that  the  basic  ideas  here  were  much 
broader  than  just  their  applications  to  microscopy.  For  the  first  time,  we 


could  produce  a  slm,  'e  optical  range  sensor  for  robotics.  We  could  also  use 
these  techniques  to  carry  out  nondestructive  measurements  of  machined  parts 
during  and  after  manufacture  and  to  determine  the  presence  of  Imperfections. 

In  the  next  two  sections,  we  will  describe  some  of  these  results  and  point  out 
their  Importance  to  nondestructive  testing  and  to  manufacturing  technology. 


Range  Sensor 

We  are  developing  two  types  of  type  II  scanning  optical  microscopes.  The 
first  mechanically  scans  the  object  or  the  lens,  the  second  uses  electronic 
scanning  and  direct  phase  and  amplitude  measurements  of  the  optical  signal. 

The  purpose  of  the  first  system  Is  to  develop  and  demonstrate  a  number  of 
new  types  of  measurement  techniques  which  should  enable  us  to  measure  range, 
film  thicknesses  of  films  a  few  wavelengths  thick,  profiles  of  films  and 
structures,  surface  roughness  and  slope,  and  the  curvature  of  a  surface.  We 
bellve  that,  at  a  later  stage,  these  techniques  will  enable  us  to  measure  the 
presence  and  size  of  dust  particles,  internal  material  properties  of  trans¬ 
parent  materials,  and  the  presence  md  size  of  surface  and  near-surface  de¬ 
fects,  such  as  cracks  and  dislocations. 

The  system  we  have  set  up  Is  Illustrated  In  Fig.  1.  A  laser  Illuminates 
a  microscope  objective  lens  with  a  collimated  beam.  The  collimated  beam  is 
obtained  by  expanding  the  laser  beam  In  a  telescope  and  passing  It  through  a 
hole  In  a  diaphragm.  The  beam  Is  focused  to  a  point  by  the  microscope  objec¬ 
tive  lens  and  an  image  of  the  point  illuminated  on  the  object  Is  imaged 
through  a  further  microscope  objective  lens  onto  a  pinhole.  A  detector  Is 
placed  behind  the  pinhole.  The  object  can  be  scanned  in  all  directions,  and 
scans  of  the  output  can  be  passed  Into  a  computer.  At  the  present  time,  line 
scans  In  any  direction  can  be  obtained.  A  Bragg  cell  can  be  placed  in  the 
path  to  amplitude  modulate  or  deflect  the  beam,  as  Illustrated  In  Fig.  1. 
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Fig.  1.  Mechanically-scanned  optical  microscope  used  for  optical  depth 
profiling,  or  V(z)  . 


In  this  system,  we  can  Insert  various  electronically-controlled  compon¬ 
ents.  At  the  present  time,  we  have  used  a  Bragg  cell  to  give  partial  elec¬ 
tronic  scanning  and  for  phase  contrast  Imaging.  We  have  employed  a  new  type 


of  electronical ly-controllable  phase  contrast  Zernlke  system.  In  addition,  we 
have  moved  the  object  up  and  down  periodically  with  a  ceramic  transducer  to 
give  an  rf  output  at  the  detector,  and  we  have  tried  out  the  concepts  with  a 
flexible,  electronical ly-controllable  mirror  to  obtain  the  position  of  the 
focus  with  extreme  accuracy. 


Optical  Range  Sensor 

The  basic  principle  of  the  type  II  microscope  is  shown  In  Fig.  2.  Light 
is  passed  through  a  pinhole  through  the  objective  lens  to  the  reflecting  ob¬ 
ject.  An  Image  of  the  pinhole  Is  obtained  on  the  reflecting  object;  this 
Image  Is  In  turn  Imaged  back  on  the  pinhole  and  a  detector  Is  placed  behind 
the  pinhole.  If  the  reflecting  plane  is  at  the  focus,  a  well -focused  image  Is 
obtained  at  the  pinhole  and  a  strong  signal  Is  obtained  at  the  detector.  If 
the  reflecting  plane  Is  displaced,  a  defocused  Image  Is  obtained  at  the  pin¬ 
hole,  and  relatively  little  light  passes  through  It.  Thus,  a  strong  maximum 
output  is  obtained  at  the  detector  when  the  beam  Is  focused  on  the  reflecting 
plane. 

The  system  Is  therefore  capable  of  very  good  range  resolution  and  does 
not  need  a  sharp  edge  to  obtain  the  range  resolution  as  does  a  self-focusing 
camera.  This  system  is  an  analog  to  the  scanning  acoustic  microscope  which 
uses  a  coherent  detector  In  the  back  focal  plane  of  the  lens.  We  have  de¬ 
veloped  a  nonparaxlal  theory  of  this  measurement  and  have  shown  that  the  range 
resolution  of  the  system  Is 

0.45  x 

d  (3  dB)  =  - 

1  -  cos  9  q 
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OPTICAL  RANGE  SENSORS 


REFLECTOR 


2.  A  diagram  illustrating  the  range  definition  of  a  type  II  microscope 


where  9q  Is  the  half-angle  of  the  aperture  of  the  lens.  The  nonparaxlal 
theory  which  we  have  developed  Implies  3  dB  range  resolutions  of  the  order 
of  .7  um  with  a  5100  A  argon  beam.  Results  obtained  with  this  system  are 
shown  In  Fig.  3.  The  comparison  between  theory  and  experiment  Is  excellent. 
The  theory  which  we  have  developed  Is  given  In  the  papers  In  the  appendix. 
Recently,  we  have  begun  to  develop  the  theory  for  polarized  optical  waves,  and 
we  are  beginning  to  take  account  of  aberrations  In  thin  films. 

The  advantage  of  this  type  of  range  sensor  Is  that  It  can  be  used  over  a 
fairly  large  range  by  placing  the  system  In  a  feedback  circuit  so  that  the 
lens  can  be  moved  mechanically  to  follow  the  surface  of  the  object.  This  re¬ 
quires  detecting  where  the  maximum  amplitude  occurs,  but  It  would  be  better 
and  more  accurate  to  look  for  a  zero.  We  have  done  this  In  a  demonstration 
experiment  by  vibrating  the  object  back  and  forth  at  a  frequency  n  with  a 
piezoelectric  pusher.  An  rf  signal  at  frequency  n  Is  obtained  at  the  de¬ 
tector.  This  signal  Is  zero  when  the  dc  output  Is  maximum,  l.e.,  we  obtain  a 
differential  of  the  output  signal.  We  have  developed  the  theory  to  show  that 
with  a  wide  aperture  lens,  we  should  be  able  to  measure  the  position  of  the 
surface  to  an  accuracy  of  .001  A  with  a  measurement  time  of  the  order  of 
1  msec  .  We  appear  to  be  within  a  factor  of  ten  of  this  accuracy,  and  are 
encouraged  by  the  early  results.  An  Illustration  of  this  type  of  differential 
measurement  Is  given  In  Fig.  4.  There  Is  a  sharp  minimum  In  the  ac  output,  as 
shown  In  the  lower  figure.  We  are  also  beginning  to  develop  the  technique  to 
measure  the  thickness  of  thin  films.  In  this  case,  we  observe  a  maximum  when 
the  focus  Is  placed  on  the  top  surface  and  the  lower  surface  of  a  transparent 
thin  film.  Early  results  on  a  ZnO  film,  1.6  um  thick,  confirm  the  basic 
Ideas.  We  are  now  developing  the  theory  to  account  for  the  aberrations  of  the 
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Plots  of  the  output  signal  V(z)  as  a  function  of  z  ,  the  distance 
from  the  focus  of  a  reflecting  plane  mirror  In  a  type  II  scanning 
microscope.  The  top  figure  has  a  numerical  aperture  of  0.9  ;  the 
lower  figure  has  a  numerical  aperture  of  0.66  . 


optical  system  when  a  wide-angle  beam  passes  Into  a  thin  film.  So  far,  the 
results,  theory,  and  experiments  appear  to  be  within  a  few  percent.  We  are 
further  working  on  the  theory  to  take  account  of  aberrations  and  to  eliminate 
errors . 

The  same  basic  Idea  applies  to  range  sensing  at  large  distances.  On  an¬ 
other  project,  not  supported  on  this  contract,  we  have  tried  out  these  Ideas 
for  measuring  range  at  distances  of  the  order  of  20  cm  .  They  appear  to  work 
well.  One  of  the  questions  that  arises  there  Is  whether  the  same  techniques 
can  be  used  on  rough  surfaces,  such  as  machined  surfaces.  It  was  convenient 
to  test  the  principles  In  the  microscope.  We  did  this  with  a  deliberately 
roughened  surface.  Fig.  5  shows  the  V(z)  curve  from  a  surface  roughened 
with  5  urn  grit.  The  surface  was  tilted  at  an  angle  of  5°  In  a  plane 
normal  to  the  z-axls  of  the  system.  This  tilt  was  sufficient  to  assure  that 
no  directly  reflected  wave  re-entered  the  optical  system. 

Clearly,  the  V(z)  behavior  In  specular  reflection  Is  Identical  to  the 
coherent  reflection  from  the  mirror,  although  the  amplitude  Is  66  dB  down. 
This  curve  also  Illustrates  the  very  high  signal -to-nolse  ratio  of  this  tech¬ 
nique,  approximately  90  dB  . 

The  results  Illustrate  that  the  basic  technique  Is  applicable  to  rough 
surfaces  because  we  are  Imaging  only  one  point  on  a  surface  at  a  time.  We 
have  worked  out  the  theory  on  this  subject  which  predicts  the  decrease  In  am¬ 
plitude  from  a  perfectly  random  surface.  Our  results  appear  to  be  In  fairly 
good  agreement  with  this  theory. 

We  would  prefer,  of  course,  to  vary  the  focal  position  of  the  beam  on  an 
rf  basis  rather  than  use  a  vibrating  object.  We  are  presently  examining  two 
techniques  for  this  purpose.  The  first  Involves  a  variable  Zernlke  phase 
plate  placed  In  the  optical  path  (Fig.  6).  The  system  Is  different  from  a 


conventional  phase  contrast  microscope  because  It  gives  an  rf  signal  output 
which  Is  directly  proportional  to  phase  change  with  a  sign  that  reverses  with 
phase,  and  Is  zero  for  zero  phase  change.  We  have  also  shown  that  It  has  much 
better  transverse  resolution.  The  phase  plate  Is  constructed  on  transparent 
PLZT  ceramic.  When  a  voltage  Is  applied  to  the  transparent  electrodes  on  the 
plate,  the  phase  of  the  optical  beam  passing  through  It  Is  changed.  If  the 
plate  Is  divided  Into  two  parts,  the  phase  of  the  central  portion  of  the  beam 
can  be  varied  relative  to  the  outer  portion.  If  the  phase  plate  Itself  Is 
placed  In  the  back  focal  plane  of  the  lens,  as  Illustrated  In  Fig.  6,  the  out¬ 
put  of  the  detector  Is  of  the  form 

cos($Q  +  A  cos  nt) 

where  A  Is  proportional  to  the  amplitude  of  the  voltage  applied  to  the  cen¬ 
ter  electrode  at  a  frequency  n  and  Is  proportional  to  the  displacement 
of  the  reflecting  plane  z  from  the  focus.  If  A  Is  small,  the  signal  obtained 
at  the  output  of  the  detector  has  a  term  that  varies  as  A  cos  nt  sin  . 

Thus  an  ac  signal  Is  obtained  which  Is  zero  when  the  focus  Is  at  the  surface 
of  the  reflecting  object.  We  have  made  crude  demonstrations  of  an  early  ver¬ 
sion  of  this  device;  It  works.  We  are  now  constructing  a  version  of  it  which 
Is  more  precisely  made  to  make  measurements  with  better  sensitivity  and  ac¬ 
curacy. 

Spherical  Mirror 

A  second  system  has  been  constructed  using  a  flexible  spherical  mirror  In 
the  path;  It  Is  Illustrated  In  Fig.  7.  This  mirror  Is  composed  of  two  pieces 
of  PZT  ceramic  bonded  together  with  a  microscope  cover  slip  on  which  a  metal 
film  Is  deposited  and  bonded  to  them.  When  voltages  are  applied,  as  shown  in 
the  figure,  one  PZT  disk  expands  and  the  other  contracts,  thus  causing  the 
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Fig.  7.  A  schematic  of  the  adjustable  mirror  system 


mirror  to  adopt  a  spherical  shape.  This  moves  the  focus  In  and  out  so  that  If 
rf  signals  are  applied  to  the  mirror  configuration,  an  rf  output  Is  obtained 
In  much  the  same  manner  as  the  phase  contrast  Zernlke  system. 

We  have  tried  this  system  out  on  our  large-scale  range  sensor  with  excel¬ 
lent  results.  In  that  situation  we  were  able  to  obtain  a  sensitivity  of  the 
order  of  50  um  at  a  distance  of  20  cm  .  This  definition  should  Improve 
with  time.  Our  theory  Indicates  that  the  definitions  In  the  microscope,  using 
either  this  system  or  the  Zernlke  phase  contrast  system,  should  be  of  the  or¬ 
der  of  0.001  A  . 

Transverse  Resolution 

The  transverse  resolution  of  these  kinds  of  systems  Is  also  very  good. 
This  Is  because  the  lens  Is  used  twice,  both  for  transmission  and  reception, 
so  that  the  response  of  the  system  Is  the  square  of  that  of  the  conventional 
microscope.  This  Implies  that  sldelobe  levels  are  low  and  the  3  dB  defini¬ 
tion  In  the  transverse  direction  Is  Improved  by  about  a  factor  of  1.4  .  We 
will  deal  with  some  of  these  basic  results  when  we  describe  the  electronic- 
ally-scanned  microscope. 

We  are  also  concerned  with  the  use  of  these  techniques  to  measure  slope 

and  curvature.  When  the  Bragg  cell  In  the  system  Is  driven  with  a  signal  of 
jugt 

the  form  e  cos  at  ,  1.e«,  with  two  signals  at  frequencies  ug  -  a  and 
u»g  +  a  ,  respectively.  It  generates  two  closely-spaced  spots  on  the  sample 
(Figs.  1  and  8).  The  signals  corresponding  to  the  two  spots  pass  along  the 
same  path,  through  the  pinhole,  to  the  detector.  If  the  spots  strike  surfaces 
with  different  heights,  they  will  generate  a  phase  signal  of  the  form 
co$(4at  +  <*)  ,  where  $  *  2kz  Is  the  phase  difference  corresponding  to  a 
displacement  z  of  one  surface  with  respect  to  the  other  and  k  =  2*/x  where 
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x  Is  the  wavelength  of  the  optical  beam.  Thus,  the  output  obtained  Is  the 
derivative  of  the  phase  variation  along  the  sample,  l.e..  It  Is  the  slope  of 
the  surface.  If  the  beams  pass  over  a  sharp  step,  a  large  peak  In  output  Is 
obtained  when  they  are  equally  spaced  from  the  step.  This  makes  It  possible 
to  locate  an  edge  with  great  accuracy.  An  Illustration  of  such  results  Is 
shown  In  Fig.  9. 

We  suggest  that,  by  using  a  modulated  signal  with  three  frequencies  In 
It,  we  can  produce  three  spots  and  therefore  use  the  system  to  measure  curva¬ 
ture  rather  than  slope.  This  could  be  extremely  useful  for  measuring  the 
diameter  of  a  hole,  the  eccentricity  of  a  round  object,  or  the  wear  of  an 
engine  cylinder. 

We  are  trying  yet  other  techniques  which  Involve  measurements  of  ac  flue 
tuatlons  associated  with  photoacoustic  modulation  by  a  modulated  laser  beam 
using  these  methods.  They  are  very  powerful  for  measuring  thermal  bonds  and 
Internal  defects, 

ELECTRONICALLY-SCANNED  OPTICAL  MICROSCOPE 

During  the  last  six  months  we  have  made  preliminary  measurements  on  the 
system  built  Into  a  Leltz  microscope.  The  Leltz  Company  has  loaned  us,  on 
Indefinite  loan,  a  modified  orthoplan  microscope  with  provision  made  to  Image 
the  back  focal  plane  of  the  objective  lens  on  a  plane  outside  the  microscope 
where  the  Bragg  cell  can  be  mounted.  The  light  from  an  argon  laser  Is  passed 
through  a  fiber-optic  delay  line  to  the  Bragg  cell. 

The  basic  system  Is  Illustrated  In  Fig.  10.  Transfer  lenses  are  mounted 
inside  the  microscope,  as  shown  In  Fig.  11.  At  the  same  time.  It  Is  possible 
to  use  the  microscope  In  the  normal  manner  with  a  separate  light  source  and 
look  at  the  object  being  measured.  This  Is  a  major  Improvement  because  It 


Fig.  10.  Basic  configuration  of  the  electronically-scanned  microscope 


SCANNED 

BEAM 


REFERENCE 

BEAM 


Fig.  11.  Schematic  of  electronically-scanned  microscope  with  transfer  lenses 
and  a  field  step.  The  second  set  of  beams  (tilted  Into  the  paper)  Is 
not  shown. 


makes  the  problems  of  lining  up  the  system  and  deciding  on  a  region  of  In¬ 
terest  on  a  sample  far  simpler. 


Fast  Oata  Acquisition 


We  have  recently  completed  the  construction  and  testing  of  an  advanced 
data  acquisition  system.  Completely  computer  controlled*  It  Includes  a  dig¬ 
itally  programmable  VCO  to  provide  scanning  on  a  random  access  basis  (we  can 
access  any  points  on  a  line  In  any  order  we  want  to),  a  logarithmic  amplitude 
detector  with  12-bit  accuracy,  80  dB  dynamic  range,  and  most  prominently,  a 


successive-approximation  phase  digitizer  with  an  accuracy  of  1/10  degree 
(corresponding  to  an  uncertainty  In  sample  height  of  about  one  Angstrom).  All 
three  are  capable  of  operating  together  at  a  rate  of  50,000  to  70,000 
points  per  second. 


This  represents  an  Improvement  of  500  times  In  speed  and  3-10  times 
In  accuracy  over  our  previous  system,  which  was  based  on  a  lock-in  amplifier. 
With  Improved  computer  hardware  and  software.  In  particular  DMA  transfer  of 


instructions  and  data,  we  could  take  512  point  lines  at  100-140  lines  per 
second . 


Submicron  Lines 

We  have  made  preliminary  measurements  with  the  Leltz  microscope  while 
building  the  fast  data  acquisition  system.  The  worst  problem  that  we  have 
encountered  with  the  present  system  Is  a  poor  signal -to-nolse  ratio.  This  Is 
mainly  due  to  loss  of  light  In  the  beam  splitters  Inside  the  microscope.  The 
one-»ay  loss  of  each  beam  Is  of  the  order  of  7  dBs  ,  but  since  the  beams  are 
used  In  a  reflection  mode  and  their  product  Is  taken  In  the  detector,  the  out¬ 
put  power  from  the  detector  Is  decreased  by  28  dB  ,  a  serious  loss.  This 
will  be  remedied  by  replacement  of  one  of  the  beam  splitters,  the  main  source 


of  loss,  and  by  Improvement  In  the  broadband  detector  matching  circuit.  We 
expect  to  obtain  a  30  dB  Improvement  In  slgnal-to-nolse  ratio  with  these 
stratagems  • 

We  have  carried  out  scans  of  900  A  thick  step  of  aluminum  on  alum* 

Inum.  The  measured  phase  variation  across  this  step  Is  shown  In  Fig.  12.  A 
theoretical  curve  is  plotted  on  top  of  the  experimental  results.  The  agree¬ 
ment  between  experiment  and  theory  Is  relatively  good.  In  these  results,  the 
distance  between  the  10-90%  height  points  Is  approximately  0.25  um  (0.5 
wavelengths);  thus,  the  resolution  of  the  system  Is  already  very  good,  and 
better  than  with  the  standard  optical  microscope. 

Because  we  have  both  amplitude  and  phase  Information  available,  we  can 
Inverse  filter  the  data  to  obtain  better  transverse  resolution.  When  we  do 
this  using  different  choices  for  the  deconvolution  filter,  we  obtain  the  re¬ 
sults  shown  In  Fig.  12.  The  choice  of  deconvolution  filter  Is  dictated  by  how 
much  ripple  we  wish  to  obtain  at  the  edges  of  the  step:  the  worse  the  ripple, 
the  worse  the  resolution.  An  Intermediate  curve  shown  gives  a  10-90%  edge 
definition  of  the  order  of  0.15  um  after  post-processing.  We  have  carried 
out  similar  measurements  on  a  1000  A  thick  metal  film  0.97  um  (1.8 
wavelengths)  wide.  Referring  to  Fig.  13,  the  upper  figure  shows  the  uncor¬ 
rected  measurement  of  the  film;  we  observe  both  Its  height  and  its  width. 

After  Inverse  filtering,  the  shape  becomes  almost  rectangular,  and  the  height 
of  the  curve  Is  almost  at  Its  full  value  for  a  step.  Therefore,  this  proced¬ 
ure  enables  us  to  measure  simultaneously  both  height  and  width  of  the  film,  a 
major  Improvement  over  most  current  metrology  systems. 

The  high  performance  Is  further  Illustrated  by  the  sharp  Image  of  the 
0.61  pm  (1.2  wavelengths)  wide  line,  shown  In  Fig.  14.  Micron  and  submicron 
line  widths  are  generally  measured  by  fitting  a  model  to  the  observed  data  and 
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Fig.  12.  A  comparison  of  three  results  taken  on  a  900  A  step  of  aluminum  on 
aluminum  with  the  electronically-scanned  microscope.  The  experiment¬ 
al  results  are  compared  with  theory,  along  with  improvements  In  the 
step  response  made  by  using  two  different  Inverse  filters. 
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Fig.  13.  A  comparison  of  phase  measurements  taken  on  a  1000  A  thick  metal 
film  0.97  urn  wide  (upper  figure),  with  the  results  (lower  figure) 
after  Inverse  filtering. 
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extracting  height  and  width  from  the  model  parameters.  The  great  sharpness 
and  the  amplltude/phase  capabilities  of  our  system  make  possible  much  more 
accurate  measurements  by  presenting  more  Information  to  the  model. 


Further  Developments 

In  the  next  few  weeks,  we  plan  to  make  two  improvements  In  the  optics  of 
the  system.  The  present  system  has  four  beams  Incident  on  the  sample,  which 
necessitates  the  use  of  special  samples  containing  a  featureless,  flat,  mir¬ 
rored  area  for  two  of  the  beams.  The  second  set  of  beams  Is  used  to  produce  a 
reference  electronic  signal.  We  plan  to  reflect  this  second  pair  of  reference 
beams  from  a  removable  plane  mirror  located  at  an  Internal  field  stop  of  the 
microscope,  located  at  an  Image  plane  of  the  object,  as  shown  In  Fig.  11. 

This  mirror  duplicates  the  function  of  the  special  samples  without  their  dif¬ 
ficulties. 

In  addition,  we  plan  to  Introduce  a  second,  fixed  acoustic  frequency  Into 
the  Bragg  cell.  This  will  produce  a  second  stationary  spot  near  one  edge  of 
the  scanned  region  which  will  take  over  the  function  of  the  undiffracted  spot 
In  the  present  system  without  the  attendant  loss  of  fleld-of-vlew.  In  ad¬ 
dition,  the  extra  frequency  shift  will  greatly  simplify  the  rf  signal  proces- 
i  sing. 

CONCLUSION 

We  have  demonstrated  a  number  of  new  measurement  systems  suitable  for 
long-distance  range  sensing  as  well  as  microscopic  range  and  transverse  resol¬ 
ution  sensing.  The  systems  are  applicable  for  use  with  photolithography,  for 
material  measurements,  and  for  measurements  during  manufacture,  as  well  as  In 
I  applications  In  nondestructive  testing  of  small  and  large  samples  In  the 

field. 
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Abstract— Tttt  measurement  and  the  use  of  phase  in  acoustic  micros¬ 
copy  are  discussed.  It  is  demonstrated  that  in  many  applications  phase 
can  be  used  to  provide  sensitivity  and  information  unparalleled  by  am¬ 
plitude-only  measurement  methods.  A  technique  capable  of  high- 
accuracy  measurement  of  the  phase  of  short  RF  acoustic  pulses  is  de¬ 
scribed.  The  power  of  this  phase  measurement  technique  is  illustrated 
in  a  number  of  applications.  Surface  material  property  measurements 
such  as  the  Rayleigh-wave  velocity  and  the  inversion  of  the  complex  T(;| 
to  obtain  the  reflectance  function  of  a  liquid-solid  interface  are  consid¬ 
ered.  Surface  topography  mapping  based  on  phase  measurement  is  ex¬ 
amined.  A  Fourier  transform  approach  for  precision  determination  of 
linewidths  comparable  to  the  resolution  spot  size  is  also  presented. 

I.  Introduction 

HE  SCANNING  acoustic  microscope  is  a  high-reso¬ 
lution  imaging  system  in  which,  unlike  conventional 
optical  microscopes,  it  is  relatively  straightforward  to 
measure  the  phase  of  the  return  signal.  With  only  a  few 
exceptions  1 1  )-{3J,  attention  has  centered  only  on  the  use 
of  intensity  information.  In  some  of  the  main  areas  of  ap¬ 
plication  of  the  acoustic  microscope,  the  phase  of  the  re¬ 
ceived  signal  plays  an  important  role.  For  example,  in  the 
so-called  F(r)  measurements  [4|-|7|,  the  amplitude  of  the 
received  signal  F(;)  as  a  function  of  the  separation  be¬ 
tween  the  lens  and  the  substrate  exhibits  periodic  peaks 
and  nulls.  This  phenomenon  is  due  to  the  beating  between 
a  specularly  reflected  signal  from  the  substrate  and  a  de¬ 
layed  leaky  Rayleigh-wave  signal,  which  reemits  to  the  lens 
while  propagating  along  the  surface  of  the  substrate.  The 
phase  difference  between  these  two  signal  components  is 
therefore  of  great  importance,  and  in  fact  it  controls  the 
contrast  of  rcflectional  images.  In  addition,  it  has  been 
shown  that  the  independent  measurement  of  phase  and 
amplitude  can  be  very  useful  in  the  determination  of  the 
clastic  constants  of  tissue  [8],  [9|.  However,  little  effort 
has  been  made  to  extract  the  phase  of  the  return  signal 
separately. 

One  reason  for  the  reluctance  to  make  use  of  the  phase 
is  that  it  is  not  generally  trivial  to  measure  the  phase  of  a 
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high-frequency  tone  burst  with  sufficient  accuracy.  We  will 
describe  here  an  approach  that  is  relatively  easy  to  imple¬ 
ment  and  capable  of  yielding  high-precision  phase  data, 
even  from  very  short  tone  bursts.  The  details  of  the  phase 
measurement  technique  will  be  published  elsewhere,  and 
we  will  concentrate  here  on  some  of  the  applications  of 
the  system  to  acoustic  microscopy. 

There  are  many  motivations  for  making  measurements 
of  this  type.  Firstly,  they  offer  a  direct  indication  of  ma¬ 
terial  properties  as  in  the  measurement  of  Rayleigh-wave 
velocity.  Also,  the  results  of  the  phase  measurement  can 
be  used  to  infer  the  width  and  the  height  of  surface  fea¬ 
tures.  Finally,  the  combination  of  amplitude  and  phase 
measurements  can  be  used  in  the  reconstruction  processes 
in  which  complete  information  about  the  interaction  be¬ 
tween  acoustic  field  and  material  is  required,  as  in  the  in¬ 
version  of  F(’)  data  to  find  the  reflection  coefficient  as  a 
function  of  angle  (10). 

II.  Phase  Measurement  Scheme 
A.  Acoustic  Lens  Configuration 

In  making  precise  phase  measurements  that  are  related 
to  physical  properties,  it  is  generally  important  to  ensure 
that  the  reference  signal,  against  which  the  phase  of  the 
probing  signal  is  to  be  compared,  and  the  probe  itself  share 
as  many  of  the  instrumental  and  environmental  phase  dis¬ 
turbances  as  possible.  In  other  words,  precise  measure¬ 
ments  are  best  done  in  an  interferometer,  where  the  two 
arms  are  closely  matched. 

In  this  work  we  have  made  use  of  two  acoustic  mea¬ 
surement  configurations  that  largely  satisfy  this  general 
condition.  The  first  (Fig.  1(a))  is  reminiscent  of  the  de- 
focused  condition  used  in  F(z)  measurements,  except  that 
here  a  fixed  separation  between  the  lens  and  the  sample 
surface  is  maintained.  Two  components  of  the  acoustic 
field  returned  to  the  lens  contribute  most  significantly  to 
the  output  signal:  the  on-axis  specular  reflection  of  the 
longitudinal  wave  in  the  water  (L  in  Fig.  la);  and  the  off- 
axis  rays  (R),  which  satisfy  the  condition  vwlvR  =  sin  dR 
for  the  conversion  of  longitudinal  waves  with  velocity  vw 
in  the  water  to  leaky  Rayleigh  waves  with  velocity  vR  on 
the  surface  of  the  sample  [41,  [5] . 

These  two  signals,  L  and  R ,  differ  slightly  in  path  length 
in  water  but  otherwise  experience  much  the  same  environ- 
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Fig.  I.  Acoustic  lens  configurations,  (a)  Rayleigh-wave  velocity  perturba¬ 
tion  measurement,  (b)  Topography  mapping. 

mental  disturbances.  However,  the  Rayleigh- wave  com¬ 
ponent  acquires  substantial  phase  delay  at  the  surface  of 
the  specimen,  hence  the  effective  path  length  for  the  R 
signal  is  longer  than  that  of  the  specular  reflection.  Indeed 
in  our  experiments  we  use  the  different  arrival  time  of 
these  two  return  signals  to  facilitate  their  separation  and 
subsequent  phase  comparison.  The  excitation  signal  is  two 
to  three  cycles  in  duration,  and  the  defocus  distance  is 
sufficient  that  there  is  no  temporal  overlap  of  the  signals. 
It  is  clear  that  changes  in  the  Rayleigh-wave  velocity  as  a 
function  of  position  along  the  surface  of  the  specimen  may 
be  sensed  in  this  way.  By  combining  a  high-accuracy  phase 
measurement  scheme  with  a  suitably  extended  Rayleigh- 
wave  path,  remarkably  sensitive  measurements  can  be 
made. 

The  second  configuration  that  concerns  us  is  illustrated 
in  Fig.  1(b).  In  this  case  the  lens  is  positioned  so  that  its 
focus  is  at  the  surface  of  the  specimen  or  slightly  above  it 
so  that  no  Rayleigh  wave  of  importance  is  excited.  A  ref¬ 
erence  path  is  provided  by  an  annular  beam  that  propa¬ 
gates  through  the  flat  outer  periphery  of  the  lens.  By  ex¬ 
citing  the  lens  with  a  short  pulse,  and  once  again  using 
time  discrimination  to  separate  the  signals  from  the  two 
different  paths,  an  interferometer  is  formed.  Now  small 
local  changes  in  the  surface  topography  result  in  large 
changes  in  the  phase  of  the  focused  beam  relative  to  the 
phase  of  the  large  diameter  reference  beam. 

In  both  measurement  configurations  it  is  apparent  that 
the  signals  of  interest  arrive  at  different  times,  and  this 
would  ordinarily  make  phase  comparison  impossible. 
However,  the  electronic  system  described  in  the  next  sec¬ 
tion  essentially  reconstructs  two  phase-coherent  continu¬ 
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Fig.  2.  Schematic  for  the  phase  measurement  system. 


ous  signals  that  are  directly  related  to  the  amplitudes  and 
phases  of  the  two  time-distinct  return  signals. 

The  first  configuration  is  also  used  where  the  signal  of 
interest  is  the  transducer  output  voltage,  when  the  system 
is  operated  in  the  F(z)  mode.  In  this  case  the  exciting 
signal  has  a  much  longer  duration  so  that  the  return  signals 
( L  and  R)  overlap  and  interfere  at  the  transducer  even  for 
large  defocus  distances.  The  idea  here  is  to  measure  both 
the  amplitude  and  the  phase  of  the  signal  relative  to  some 
reference  as  a  function  of  the  axial  separation  z  between 
the  acoustic  lens  and  the  specimen.  The  resulting  complex 
V(z)  can  be  inverted  1 101  to  obtain  the  angular  depen¬ 
dence  of  the  reflectance  of  the  liquid-specimen  interface. 
From  the  reflectance  function,  various  material  property 
parameters  can  readily  be  extracted. 

In  the  experiments  described  here,  a  center  frequency 
of  50  MHz  was  used  with  an  acoustic  transducer  having  a 
bandwidth  of  20  percent.  The  lens  material  was  fused 
quartz  and  had  a  radius  of  curvature  of  3.2  mm,  giving  a 
focal  length  in  water  of  4.24  mm.  The  opening  aperture 
of  the  lens  was  5.0  mm  in  diameter,  corresponding  to  a 
maximum  half  angle  of  36°  or  an  f-number  of  0.85. 

Providing  the  same  fractional  bandwidth  can  be  main¬ 
tained,  there  should  be  little  difficulty  in  applying  the  same 
techniques  to  a  system  operating  at  much  higher  fre¬ 
quency,  perhaps  up  to  2  GHz.  Beyond  that  frequency,  lim¬ 
itations  in  the  electronic-switching  components  presently 
available  may  present  difficulties. 

B.  Electronics 

The  essential  elements  of  the  phase  measurement  elec¬ 
tronics  are  shown  schematically  in  Fig.  2.  The  system  will 
be  described  in  detail  in  a  forthcoming  paper  (18],  where 
various  sources  of  error  and  an  analysis  of  the  ultimate 
performance  are  discussed.  A  more  general  review  of  the 
concepts  will  suffice  here. 
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The  basic  scheme  involves  signal  recovery  by  synchro¬ 
nous  detection  and  subsequent  phase  measurement  with  a 
lock-in  amplifier.  The  key  component  (Fig.  2)  in  the  signal 
source  is  the  single  sideband  generator  (SSB)  with  syn¬ 
chronous  outputs  at  100  kHz,  10.6  MHz,  and  10.7  MHz. 
The  10.7-MHz  output  is  in  fact  the  upper  sideband  of  the 
product  of  the  100  kHz  and  10.6-MHz  signals.  The  lower 
sideband  at  10.5  MHz  is  suppressed  by  at  least  50  dB  using 
a  standard  FM  radio  system  IF  filter.  The  translation  os¬ 
cillator  shifts  the  operating  frequency  up  to  the  desired 
center  frequency  of  the  acoustic  system,  which  is  50  MHz 
in  this  case.  A  reference  signal  at  49.9  MHz  is  also  gen¬ 
erated.  Again  the  lower  sideband  components  of  the  prod¬ 
uct  signals  have  to  be  removed.  At  this  point  however  the 
filtering  requirements  are  not  as  stringent  as  in  the  SSB 
generator  and  tunable  bandpass  filters  (F3  and  F4)  with 
five-percent  bandwidth  are  adequate,  since  the  sidebands 
are  now  widely  separated.  This  signal  generation  scheme 
is  flexible  in  that  the  operating  frequency  is  tunable  over  a 
fairly  wide  range,  which  is  limited  in  this  case  to  100  MHz 
by  the  bandpass  filters.  For  much  higher  frequency  of  op¬ 
eration,  filters  with  greater  selectivity  would  be  needed  or, 
alternatively,  successive  stages  of  heterodyning  could  be 
used  with  filtering  at  each  stage  to  step  up  incrementally 
to  the  desired  operating  frequency. 

As  illustrated  in  Fig.  2  with  a  50-MHz  acoustic  micro¬ 
scope,  the  50-MHz  continuous  wave  (CW)  signal  is  time¬ 
gated  to  produce  a  short  tone  burst  that  excites  the  acous¬ 
tic  transducer.  Switch  SW3  is  used  as  a  time  gate  to  pass 
the  low  -level  acoustic  return  signals  of  interest  and  to  block 
high  level  extraneous  ones,  which  may  damage  the  pream¬ 
plifier.  The  acoustic  return  signals  are  time-separated  RF 
pulses,  and  they  are  electronically  separated  through  time¬ 
gating  into  two  channels  as  shown.  Each  of  the  resulting 
signals  is  mixed  with  the  49.9  MHz  reference  signal,  and 
the  product  is  narrow-band  filtered  to  extract  the  100-kHz 
component. 

It  can  be  shown  118]  that  the  100-kHz  signals  are  es¬ 
sentially  low-frequency  CW  replicas  of  the  pulse  modu¬ 
lated  RF  acoustic  signals,  bearing  identical  phase  and  am¬ 
plitude  information.  With  a  biphase  lock-in  amplifier  tuned 
to  100  KHz,  one  can  readily  measure  the  phase  difference 
between  the  two  signals  and  also  the  amplitude  of  the  sig¬ 
nal  being  fed  to  the  "signal”  channel  of  the  lock-in,  if  it 
is  required.  Phase  sensitivity  is  basically  limited  by  the 
lock-in  amplifier  because  the  problem  of  system  noise  can 
in  most  cases  be  overcome  by  increasing  the  integration 
time  in  the  lock-in  output  stage.  This  effectively  reduces 
the  system  noise  bandwidth,  but  it  also  includes  the  ob¬ 
vious  disadvantage  of  longer  data  acquisition  time.  Good- 
quality  lock-in  amplifiers  routinely  have  phase  resolution 
of  the  order  of  0.1°,  which  corresponds  to  an  overall  sys¬ 
tem  phase  sensitivity  of  1/3600  of  a  wavelength. 

In  velocity  perturbation  measurements  the  measured 
phase  is  dependent  on  the  distance  between  the  lens  and 
the  specimen.  This  distance  changes  with  the  sample  sur¬ 
face  topography  and  is  susceptible  to  thermal  drift  as  well. 


thus  introducing  phase  error  into  the  measurement.  To 
minimize  this  error,  a  feedback  mechanism  is  employed 
to  keep  the  lens-to-specimen  spacing  constant.  The  acous¬ 
tic  lens  is  mounted  on  a  piezoelectric  (PZT)  stack  so  that 
its  vertical  position  can  be  adjusted  continuously  by  an 
electronic  control  signal.  The  acoustic  on-axis  reflection 
pulse,  whose  phase  is  a  direct  measure  of  the  lens  to  sam¬ 
ple  distance,  is  applied  to  the  reference  channel  of  the  lock- 
in  amplifier  as  shown  in  Fig.  2.  The  reference  channel 
generates  a  100-kHz  constant-amplitude  phase-locked 
replica  of  the  reference  input,  which  is  then  compared  with 
the  100-kHz  output  of  the  SSB  generator  to  produce  the 
control  signal  for  the  PZT  stack.  The  acoustic  lens  auto¬ 
matically  tracks  the  surface  topography  of  the  sample  dur¬ 
ing  scanning  to  ensure  that  the  measured  phase  change  is 
due  to  material  property  variation  alone.  This  feedback 
mechanism  is  also  used  effectively  in  such  modes  of  op¬ 
eration  as  topography  measurement  to  compensate  for 
thermal  effects. 

III.  Applications 

The  power  of  this  measurement  scheme  is  illustrated  in 
the  following  with  a  number  of  examples.  We  will  consider 
both  material  property  measurements,  such  as  Rayleigh- 
wave  velocity  and  complex  reflectivity,  end  surface  topog¬ 
raphy  measurements. 

A.  Velocity  Perturbation  Measurements 

In  this  section  we  will  specifically  deal  with  the  pertur¬ 
bation  to  the  Rayleigh-wave  propagation  velocity  as  a  re¬ 
sult  of  material  property  change  and  also  the  presence  of 
surface  residual  stress. 

I)  Velocity  Perturbation  Due  to  Thin- Film  Over¬ 
lay:  The  sample  in  this  example  is  a  multiple-thickness 
indium  film  deposited  on  glass.  The  thicknesses  are  240 
A  and  620  A ,  respectively.  The  objective  was  to  measure 
the  perturbation  of  the  Rayleigh-wave  velocity  caused  by 
the  indium  film.  A  line  scan  over  the  surface  of  the  sample 
(Fig.  3)  exhibits  phase  changes  of  7°  and  11 0  for  the  240- 
A  and  380- A  step  changes  in  film  thickness.  The  spatial 
resolution  of  the  system  is  defined  by  the  Rayleigh-wave 
path  length  on  the  substrate  and  is  determined  to  be  about 
0.8  mm  from  the  step  transition  width  in  the  line  scan.  It 
can  be  calculated  from  first  order  perturbation  theory  1 1 1  ] 
that  the  velocity  perturbation  due  to  the  indium  film  is  0.18 
percent  and  0.46  percent  for  the  240- A  and  620- A  layers, 
respectively.  Based  on  these  estimated  parameters,  one 
would  expect  phase  changes  of  7.6°  and  12°  for  the  240- 
A  and  380- A  step  transitions.  Hence  there  is  fairly  good 
agreement  between  the  experimentally  obtained  and  the¬ 
oretically  predicted  phase  changes.  It  should  also  be  noted 
that  in  the  line  scan,  the  small  phase  variations  in  the  sup¬ 
posedly  flat  regions  of  the  indium  film  are  real  and  re¬ 
peatable.  The  fluctuations  are  less  than  0.5°  and  are  due 
to  nonuniformity  in  the  thickness  of  the  indium  film.  Since 
the  phase  sensitivity  of  the  system  is  limited  by  the  lock- 
in  amplifier  to  0.1°,  this  measurement  technique  can  po- 
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Fig.  3.  Measurement  of  velocity  perturbation  due  to  indium  film  on  glass. 

tentially  detect  velocity  perturbations  on  the  order  of  one 
part  in  I05. 

2)  Residual  Stress  Measurement:  In  nondestructive 
testing,  knowledge  of  the  surface  residual  stress  distribu¬ 
tion  in  a  component  plays  an  important  role  in  the  predic¬ 
tion  of  failure  modes.  One  way  of  characterizing  residual 
stress  is  by  measuring  the  acoustic  wave  propagation  ve¬ 
locity,  which  varies  linearly  with  the  local  stress  [12].  With 
the  measurement  configuration  shown  in  Fig.  1(a)  the 
measured  change  in  the  relative  phase  between  the  L  and 
R  signal  pulses  as  the  lens  is  scanned  can  be  shown  to  be 
directly  proportional  to  the  residual  stress  on  the  object 
surface  [13]. 

The  sample  used  in  the  experiment  is  a  glass  disk  of  two 
inches  in  diameter  that  was  heated  and  then  thermally 
quenched  with  air  jets  to  introduce  a  radial  distribution  of 
residual  stress.  Fig.  4(b)  shows  the  radial  variation  of  the 
measured  phase  perturbation.  For  comparison  a  destruc¬ 
tive  test  was  carried  out  on  a  similar  sample  to  determine 
the  actual  residual  stress  distribution.  A  Vicker’s  indenter 
was  used  at  a  prescribed  load  to  produce  median  cracks  at 
different  points  on  the  glass  disk.  By  measuring  the  crack 
lengths,  the  residual  stress  as  a  function  of  radial  distance 
was  calculated  [14].  The  result  of  the  destructive  test  is 
shown  in  Fig.  4(a).  The  variation  of  the  principal  stresses 
are  essentially  the  same,  with  the  radial  and  tangential 
components  tracking  each  other  to  within  20  MPa.  Since 
the  phase  measurement  is  omnidirectional  in  that  a  spher¬ 
ical  lens  was  used  to  launch  Rayleigh  waves  propagating 
in  all  directions  along  the  surface,  the  resulting  phase  per¬ 
turbation  is  a  measure  of  the  sum  of  the  principal  stresses. 
Comparison  of  Figs.  4(a)  and  4(b)  shows  that  there  is 
fairly  good  corroboration  between  the  phase  perturbation 
curve  and  the  actual  residual  stress  distribution.  From 
these  results  we  empirically  deduce  that  40  MPa  of  stress, 
which  is  the  change  from  the  center  of  the  sample  to  the 
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Fig.  4.  Residual  stress  measurements  on  the  surface  of  a  tempered  glass 
disk,  (a)  Variation  of  radial  (open  symbols)  and  tangential  stresses  (closed 
symbols)  across  a  diameter  of  the  glass  disk  estimated  from  indentation 
fracture  tests,  (b)  Radial  variation  of  measured  phase  perturbation. 

limit  of  the  line  scan,  corresponds  to  a  phase  variation  of 
30°. 

B.  Surface  Topography  Measurements 

We  have  used  the  acoustic  microscope  in  the  measure¬ 
ment  configuration  shown  in  Fig.  1(b)  as  a  high-resolution 
noncontacting  profilometer  [13],  [15],  Besides  being  able 
to  make  use  of  phase  profile  to  map  the  depth  variation, 
we  can  also  determine  with  high  precision  the  transverse 
profile  of  surface  features  in  special  cases.  Of  particular 
interest  in  this  regard  is  the  measurement  of  the  widths  of 
long  rectangular  strips,  which  is  a  matter  of  great  concern 
in  the  fabrication  of  semiconductor  components.  We  will 
show  here  both  experimentally  and  theoretically  the  ad¬ 
vantages  of  utilizing  phase  to  determine  the  linewidth.  We 
will  also  introduce  here  a  Fourier  transform  technique  for 
linewidth  measurement  in  cases  where  the  strip  width  is 
comparable  to  the  spot  size. 

1)  Depth  Profiling:  Topography  images  of  metallized 
stripe  patterns  on  a  fused-silica  substrate  are  shown  in 
Figs.  5(a),  5(b)  and  5(d).  The  patterns  have  progressively 
finer  pitches  and  the  linewidths  are  250,  125,  and  62.5 
fun,  respectively.  The  metallization  is  gold  with  a  thick¬ 
ness  of  about  3000  A .  The  gold  stripes  show  up  as  bright 
areas  in  the  images.  At  50  MHz,  with  an  f-number  of  0.85 
and  a  uniformly  excited  aperture,  the  Rayleigh  resolution 
of  the  acoustic  lens  is  1.13  FX  =  0.96  X  or  29  /tm,  while 
the  3-dB  resolution  is  0.64  FX  =  0.54  X  or  16  nm.  The 
62.5-^m  line  image  in  Fig.  5(d)  is  clearly  resolved,  as 
would  be  expected.  The  amplitude  image  of  the  \25-pm 
stripes  is  shown  in  Fig.  5(c)  for  comparison.  There  is  vir- 
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Fig.  5.  Topography  images  of  striped  metallization  patterns  on  a  fused 
artz  substrate.  The  metallization  is  gold  and  the  thickness  is  about  3000 
.  (a)  Phase  image  of  250-(im  lines,  (b)  Phase  image  of  125-pm  lines. 
(c)  Amplitude  image  of  125-pm  lines,  (d)  Phase  image  of  62.5-pm  lines, 
(e)  Perspective  plot  of  the  measured  topography  map  of  the  250-#im  line 
pattern. 

tually  no  contrast  to  suggest  the  presence  of  a  striped  pat¬ 
tern,  clearly  illustrating  that  phase  is  a  far  more  sensitive 
means  of  gauging  distances.  The  measured  phase  differ¬ 
ence  between  the  surface  of  the  metallization  and  the 
surface  of  the  fused-silica  substrate  is8.3°,  which  corre¬ 
sponds  to  a  film  thickness  of  2800  A  after  taking  into 
account  the  angular  dependence  of  the  complex  reflec¬ 
tance  function  and  the  effect  of  focusing. 

We  will  show  here  that,  based  on  theoretical  consider¬ 
ations,  indeed  the  phase  rather  than  the  magnitude  of  the 
acoustic  microscope  output  V(x)  is  a  more  sensitive  mea¬ 
sure  of  the  depth  profile.  Consider  the  simple  case  of  a 
rectangular  strip  of  width  w  on  a  substrate  of  identical 
material.  Assume  that  the  reflection  coefficient  T0  of  the 
surface  is  real.  Then  the  spatial  response  of  the  strip 
structure  is  given  by 

R(x)  =  r0  +  r0  (e»  -  1)  rect  (jc/w)  (1) 


where  0  is  the  phase  change  due  to  the  thickness  of  the 
strip.  Let  the  response  of  the  imaging  system  to  a  line 
source  be  h(x).  Furthermore,  suppose  that  the  lens  is  uni¬ 
formly  illuminated  in  the  back  focal  plane  and  the  lens  is 
aberration-free  so  that  h(x)  is  real.  The  acoustic  micro¬ 
scope  output  is  therefore  given  by 

V(x)  =  R(x)  *  h( x)  (2) 

or 

V(x)  =  r 0K  +  r0  (e1*  -  1)  s(x)  (3) 

where 

j(jt)  =  h(x)  *  rect  (x/w)  (4) 

and 

J+OO 

h(x)  dx. 

—  Oo 

The  first  term  on  the  right-hand  side  of  (3)  corresponds  to 
the  constant  background  reflection  from  the  substrate,  and 
the  second  term  corresponds  to  the  additional  spa¬ 
tially  varying  contribution  that  is  due  to  the  strip.  It  can 
easily  be  shown  that  for  a  thin  strip  where  0  «  1,  to 
second  order  in  0,  the  magnitude  and  the  phase  of  V(x) 
are  given  respectively  by  the  relations 

ymt(x)  =  r 0K  +  [s(x)  -  K )  02  (5) 

and 

^phAseW  =  £*(*)•  (6> 

The  Kphas(.(jr)  is  directly  proportional  to  0,  whereas  the 
spatially  varying  part  of  Fmag(x)  depends  on  0  to  the  sec¬ 
ond  order.  Thus  for  thin  strips,  phase  is  a  much  more  sen¬ 
sitive  measure  of  the  depth  profile. 

2)  Linewidth  Measurement:  A  major  advantage  of 
scanned  confocal  imaging  systems  such  as  ours  is  that  be¬ 
cause  the  point  spread  function  is  always  positive  and  falls 
off  rapidly  at  large  distances,  there  is  hardly  any  ripple  in 
the  amplitude  of  the  step  response  116],  117],  A  similar 
smooth  transition  in  the  phase  of  the  step  response  can 
also  be  observed  in  the  perspective  plot  shown  in  Fig.  5(e). 
For  the  purpose  of  comparison,  we  have  carried  out  a  nu¬ 
merical  analysis  for  the  step  phase  responses  of  different 
imaging  configurations.  For  small  step  heights,  or  equiv¬ 
alently  small  phase  changes,  the  normalized  phase  step 
response  is  of  the  form 

@(x)  =  |  g(r)  cos-1  (x/r)  r  dr  (7) 

where  x  is  the  distance  between  the  center  of  the  beam 
and  the  step,  and  g(r)  is  the  point  response  function  of 
the  imaging  system.  It  is  assumed  that  regions  on  either 
side  of  the  step  are  of  equal  reflectivity.  The  step  response 
function  is  plotted  in  normalized  form  in  Fig.  6  for  two 
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Fig.  6.  Step  phase  response  of  a  normal  phase  interference  microscope 
(solid  curve)  and  a  confocal  imaging  system  (theory  in  dashed  curve, 
experimental  result  in  dosed  circles). 

different  cases.  The  solid  curve  corresponds  to  g(r)  = 
jinc  (2  r  sin  0/X),  which  is  the  point  response  of  a  normal 
phase  interference  microscope  where  the  illumination  is 
a  plane  wave.  Here  6  is  the  maximum  half  angle  of  the 
objective  lens:  X  is  the  wavelength  of  the  illumination;  and 
jinc  (  r)  is  defined  as  2J^irx)/trx.  Note  that  in  this  case 
j?(  r)  reverses  in  sign  and  its  corresponding  step  phase  re¬ 
sponse  does  not  change  monotonically  with  distance.  The 
dashed  curve  corresponds  to  g(r)  =  jinc2  (2  r  sin  0/X), 
the  point  response  of  a  confocal  system.  The  confocal  re¬ 
sponse  shows  a  distinct  lack  of  ripples.  The  experimen¬ 
tally  measured  step  phase  response  is  also  plotted  in  closed 
circles  showing  close  agreement  with  theory.  In  both  im¬ 
aging  configurations,  the  50-percent  threshold  of  the  phase 
profile  demarcates  the  location  of  the  edge.  It  is  obvious, 
however,  that  the  smooth  transition  of  the  confocal  step 
phase  response  enables  a  straightforward  and  unambigu¬ 
ous  determination  of  not  only  the  strip  thickness  but  also 
the  strip  width. 

The  50-pcrcent  threshold  criterion  used  to  locate  the 
edges  of  a  strip  for  estimating  the  linewidth  is  a  valid  one 
for  wide  strips.  For  strip  widths  comparble  to  the  spot  size 
of  the  imaging  system,  the  step  responses  of  the  two  edges 
are  in  close  proximity  and  tend  to  interfere  with  each 
other,  making  the  establishment  of  a  general  criterion  for 
locating  the  strip  edges  difficult.  Wc  will  describe  here  a 
Fourier-transform  method  that  circumvents  these  difficul¬ 
ties. 

Applying  Fourier  transform  to  (6)  and  using  the  defi¬ 
nition  of  ,v(.r)  in  (4).  the  spatial  frequency  spectrum  of 
is  given  by 

*W/,)  =  |  //(/,)  *  h-  sine  ( n/()  (8) 

where  //(/,)  is  the  Fourier  transform  of  h(x).  The  term 
with  the  sine  function  corresponds  to  the  Fourier  trans- 
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Fig.  7.  Spectral  responses  of  5000- A -thick  gold  lines  on  a  fused-quartz 
substrate.  The  nominal  linewidths  arc  76  um,  43  pm.  and  33  pm.  re¬ 
spectively. 


form  of  the  spatial  variation  of  the  strip  and  its  zeroes  are 
located  at 

/,  —  n/w.  n  =  1 ,  2,  3,  ■  •  • .  (9) 

The  term  H( /,)  represents  the  angular  spectral  response 
of  the  imaging  system.  Fora  confocal  system  \  H(fx)\  is 
maximum  at/,  =0  and  generally  falls  off  monotonically. 
until  it  becomes  zero  at  the  upper  cutoff  frequency  of 

/, .  =  W(FX) 

where  F  is  the  f-number  of  the  lens,  and  X  is  the  wave¬ 
length  of  the  illumination  1 161.  1 17|.  Hence  the  zeroes  of 
Fphj,c( /,)  arc  those  of  the  strip  response  within  the  pass- 
band  of  the  imaging  system.  As  long  as 

fv  >  FX  (10) 

the  strip  width  w  can  be  determined  using  (8)  from  the 
locations  of  the  zeroes  of  FP|U.C(  /,).  In  our  case  F  =  0.85. 
so  in  principle  we  can  determine  linewidths  as  small  as 
0.85  X  or  25  pm  at  50  MHz  in  water. 

The  merit  of  this  Fourier-transform  technique  lies  in  its 
simplicity.  Given  a  priori  knowledge  of  the  strip  geome¬ 
try,  one  does  not  need  to  know  the  exact  nature  of  imaging 
system.  The  mere  location  of  the  zeroes  for  the  spectral 
response  of  the  object  suffices  for  the  accurate  determi¬ 
nation  of  the  linewidih. 

Fig.  7  illustrates  the  experimental  |Fph, „.(/,))  for  ir  = 
76  gni,  43  gm,  and  33  /oil.  respectively.  The  strips  arc 
5000- A -thick  gold  lines  on  a  fused-quartz  substrate.  The 
acoustic  microscope  operates  at  50  MHz  (X  =  30  /on) 
with  F  =  0.85,  and  the  acoustic  beam  is  focused  on  the 
surface  of  the  fused-quartz  substrate.  The  |  /, ) |  is 

plotted  on  a  log  scale.  The  weighting  due  to  the  imaging 
system  response  |  //( /) |  all  the  way  out  to  the  resolution 
limit  is  very  much  in  evidence.  The  linewidths  estimated 
from  the  position  of  the  first  zero  of  Fph  1M.< / )  are  76.3. 
43.5,  and  31.4  pm.  respectively.  They  agree  very  well  with 
the  optically  measured  widths.  The  minor  discrepancies 
are  mainly  due  to  the  nonuniformity  in  the  widths  of  the 
gold  lines. 

In  addition,  this  Fourier  transform  approach  is  essen¬ 
tially  insensitive  to  dcfocusing.  It  is  obvious  from  (8)  that 
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ITu.  H  Spoclral  responses  of  2000  A  -thick  76-pni-wide  chrome  line  on  a 
fused  quart?  substrate  al  various  deionising  distances:  ;  =  0  pm,  +60 
pm.  jnd  +  HO  pm .  respectively 


defocusing  only  changes  the  spectral  response  of  the  im¬ 
aging  system  //(/,)  slightly  and  does  not  introduce  any 
new  zeroes  to  nor  does  it  affect  the  zeroes  of  I  **«</, >1- 
Fig.  8  shows  the  experimental  |PPhase(/.)|  f°ra  2000- A  - 
thick  and  76-/cm-wide  chrome  line  on  fused  quartz  at  var¬ 
ious  defocusing  distances.  The  focus  was  located  on  the 
substrate,  and  60  /cm  and  80  /cm  above  the  substrate,  re¬ 
spectively  The  locations  of  the  first  zero  in  all  three  cases 
are  essentially  coincident,  giving  a  width  estimate  of 
76.3  /cm. 

Jt  Independent  Measurement  of  Velocity  Perturbation 
and  Topography:  In  conventional  acoustic  micrographs, 
surface  topography  and  material  property  both  contribute 
to  contrast  in  the  image,  and  their  respective  effects  are 
generally  indistinguishable  from  one  another.  With  the  to¬ 
pography  tracking  mechanism  described  in  Section  II-2 
built  into  our  phase  measurement  system,  surface  topog¬ 
raphy  and  material  property  change  can  essentially  be  ob¬ 
tained  independent  of  each  other. 

Fig.  9(a)  is  the  topographic  line  scan  of  an  aluminum 
film  step  on  glass  with  a  nominal  thickness  of  5000  A . 
The  initial  and  trailing  slow  phase  changes  are  due  to  the 
large  spatial  extent  of  the  reference  signal  as  it  traverses 
the  step,  whereas  the  abrupt  phase  change  is  due  to  the 
focused  beam  as  it  crosses  the  step.  Fig.  9(b)  is  the  same 
profile  on  an  expanded  scale.  The  phase  change  caused  by 
the  step  is  11.25°.  which  corresponds  to  a  thickness  of 
4688  A.  With  a  limiting  phase  resolution  of  0.1°,  the  ul¬ 
timate  height  sensitivity  of  the  system  is  about  50  A . 

Fig.  9(c)  shows  a  velocity  perturbation  scan  of  the  same 
sample.  The  phase  change  measured  is  7.75°,  correspond¬ 
ing  to  a  velocity  perturbation  of  0.25  percent,  which  is  in 
fair  agreement  with  the  theoretically  calculated  value  of 
0.21  percent. 

C.  Reflectance  by  Inversion  of  Complex  V(z) 

It  has  been  demonstrated  that  the  reflectance  function 
at  a  liquid-solid  interface  can  be  obtained  by  inverting  the 
corresponding  complex  F(z)  data  [10],  The  principle  mo¬ 
tivation  behind  this  endeavor  is  that  one  can  obtain  much 
more  useful  material  property  information  from  the  re¬ 
flectance  function  than  the  usual  and  rather  limited  treat¬ 
ment  of  measuring  the  periodicity  of  the  nulls  in  the  V(z) 
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Fig.  9.  Independent  measurements  of  velocity  perturbation  and  topography 
profile  due  to  a  5000- A  aluminum  film  step,  (a)  Topography  profile  scan, 
(b)  Topography  profile  scan  on  an  expanded  scale,  (c)  Velocity  pertur¬ 
bation  sen. 


curve  to  determine  only  the  Rayleigh-wave  velocity.  The 
solid  curves  in  Fig.  10  show  the  magnitude  and  the  phase 
of  the  theoretical  reflectance  function  of  a  water-fused- 
silica  interface.  The  reflectance  function  contains  infor¬ 
mation  on  the  longitudinal  and  shear  critical  angles,  and 
thus  the  velocities  of  propagation  of  the  respective  modes 
in  fused  silica.  Also,  the  Rayleigh  critical  angle  corre¬ 
sponds  to  the  point  at  which  the  phase  is  x  radians  in  the 
region,  where  the  phase  curve  undergoes  a  rapid  2x  ra¬ 
dian  change. 

It  can  be  shown  that  the  relation  between  K(z)  and  the 
reflectance  function  R(0)  is  essentially  one  of  Fourier 
transformation  [  10].  The  inversion  formula  is  given  by 

SF  {  V(u)  exp  [/'xw/sin2  (a/2)J) 


R(t)  = 


[Ut  (t)\2P2{t)[\  -  2  t  sin2  (a/2)] 


(11) 


where 


u  = 


4z  sin2  (a/2) 


and 


t  = 


sin2  (fl/2) 
sin2  (a/2)' 
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Fig.  10.  Theoretical  and  experimental  reflectance  function  for  a  water- 

fused-quartz  interface  obtained  by  inversion  of  (he  corresponding  com¬ 
plex  K(z). 

The  a  is  the  maximum  half  angle  of  the  lens;  6  is  the 
incident  angle  with  respect  to  the  liquid-solid  interface; 
and  X  is  the  wavelength  of  the  acoustic  wave  in  the  liquid 
medium.  In  addition,  Uf  is  the  back  focal  plane  illumi¬ 
nation,  and  P  is  the  pupil  function  of  the  lens.  Generally 
these  two  lens  parameters  are  not  known,  but  [t/f]2  P2  in 
the  denominator  of  (11)  can  be  calibrated  by  measuring 
the  V(z)  fora  material  whose  reflectance  function  has  uni¬ 
form  amplitude  and  phase  over  the  range  of  angular  spec¬ 
tral  components  excited  by  the  acoustic  lens. 

The  V(z)  measurements  were  conducted  at  10  MHz.  A 
water-lead  interface,  in  which  the  reflectance  function  is 
uniform  in  magnitude  and  phase  out  to  0  =  40°,  was  used 
to  calibrate  the  If/*]2  P2  term.  The  experimental  reflec¬ 
tance  function  of  fused  silica,  obtained  from  the  inversion 
procedure,  is  superposed  on  Fig.  10  in  dashed  lines.  There 
is  good  agreement  between  theory  and  experiment  in  both 
magnitude  and  phase.  The  measured  shear  critical  angle 
at  0  -  23.5°  and  the  Rayleigh  critical  angle  at  0  =  25.85° 
compare  extremely  well  with  theoretical  values. 

IV.  Conclusion 

We  have  demonstrated  that  phase  information  is  a  use¬ 
ful  asset  in  acoustic  microscopy.  Phase  can  be  used  to  pro¬ 
vide  sensitivity  and  information  in  surface  material  char¬ 
acterization  unparalleled  by  amplitude-only  measurement 
techniques.  Most  of  the  concepts  described  in  this  work 
have  complete  generality  and  can  readily  be  used  in  other 
applications.  In  fact,  many  ideas  presented  in  this  work 
have  been  applied  to  a  new  scanning  optical  microscope 
(19)  developed  by  Jungerman  and  Kino  for  surface  topog¬ 
raphy  measurement.  In  the  optical  case  both  the  amplitude 
and  the  phase  of  the  reflected  beam  can  be  measured  as 
in  the  acoustic  microscope.  Preliminary  results  from  the 
scanning  optical  microscope  give  further  confirmation  of 
most  of  the  concepts  put  foith  in  this  paper. 
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Aitimct—U  is  demonstrated  that  the  reflectance  function  R(0)  of  a 
liquid-solid  interface  can  be  obtained  by  inverting  the  complex  V(z)  data 
collected  with  an  acoustic  microscope.  The  inversion  algorithm  is  based 
on  a  non  paraxial  formulation  of  the  V(z)  integral,  which  establishes  the 
Fourier  transform  relation  between  R(0)  and  V(z).  Examples  are  given 
to  show  that  with  this  measurement  technique,  the  acoustic  phase  ve¬ 
locities  of  the  propagating  modes  in  the  solid  medium  can  easily  be 
determined  and  material  losses  can  be  estimated.  The  same  technique 
is  also  used  for  characterizing  imaging  performance  of  focused  sys¬ 
tems.  Applications  in  thin- film  measurement  are  also  discussed. 


I.  Introduction 

NEW  TECHNIQUE  for  measuring  and  interpreting 
the  acoustic  material  signature  or  the  so-called  V(z ) 
in  acoustic  microscopy  will  be  presented  in  this  paper. 
V(z)  refers  to  the  variation  of  the  normalized  voltage  out¬ 
put  of  an  acoustic  microscope  as  a  function  of  the  sepa¬ 
ration  distance  z  between  the  focal  point  of  the  acoustic 
lens  and  the  plane  surface  of  a  reflecting  object  [1]-[4J. 
This  variation  of  voltage  response  with  distance  is  char¬ 
acteristic  of  the  material  being  examined  and  is  an  im¬ 
portant  source  of  contrast  in  acoustic  micrographs.  The 
V(z)  effect  is  a  well-understood  phenomenon,  and  various 
models  have  been  put  forth  to  explain  how  it  is  related  to 
acoustic  material  properties  [1],  [2],  [5],  [6],  [9],  So  far, 
practical  application  revolves  around  the  inverse  process 
of  determining  acoustic-property-dependent  parameters  by 
interpreting  selective  features  of  the  magnitude  of  the  V(z) 
function.  For  example,  the  leaky  surface  acoustic  wave  ve¬ 
locity  can  be  obtained  directly  from  the  periodicity  of  the 
interference  peaks  and  nulls  of  the  magnitude  of  the  Viz) 
function  in  the  negative  z  region  [4]— [7],  where  the  object 
surface  is  located  between  the  focal  point  and  the  lens. 
Another  approach  involves  the  Fourier  transform  of  the 
negative  z  portion  of  the  magnitude  of  the  V{z)  function  in 
order  to  identify  propagating  modes  at  the  liquid-solid  in¬ 
terface  and  to  determine  their  respective  phase  velocities 
18). 

The  basic  theory  and  experimental  result  of  a  more  gen¬ 
eral  Fourier  transform  technique,  which  utilizes  both  the 
magnitude  and  phase  of  the  V(z)  function  for  material 
property  characterization,  will  be  given  in  this  paper.  The 
theoretical  treatment  is  a  combination  of  angular  spectra 
and  ray  optics  formalisms.  It  is  carried  out  specifically  for 
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the  case  where  the  acoustic  transducer  is  a  spherical  shell 
that  radiates  directly  into  the  liquid  coupling  medium  and 
hence  forms  a  focusing  lens  by  itself.  An  integral  expres¬ 
sion  relating  the  complex  V(z)  and  the  reflectance  R(0)  of 
the  liquid-solid  specimen  interface  as  a  function  of  the 
angle  0  between  the  direction  of  an  incident  plane  wave 
and  the  normal  to  the  interface  surface  is  derived.  It  will 
be  shown  that  by  a  suitable  transformation  of  variables, 
R(6)  can  be  obtained  by  an  inverse  Fourier  transform  of 
V(z). 

The  primary  motivation  behind  this  work  is  that  the 
complex  reflectance  function  R(0)  contains  a  wealth  of 
readily  available  acoustic  material  property  information. 
The  solid  curves  in  Fig.  7  show  the  theoretical  magnitude 
and  phase  of  the  reflectance  function  of  a  water-fused- 
silica  interface.  The  longitudinal  and  shear  critical  angles, 
6U  and  0Sc,  show  up  as  distinctive  features  on  the  magni¬ 
tude  curve  and  their  locations  give  the  propagation  veloc¬ 
ities  of  the  respective  modes  directly.  In  addition  the  Ray¬ 
leigh  critical  angle  0Kc  corresponds  to  the  point  where  the 
phase  curve  becomes  t  in  a  region  where  it  undergoes  a 
rapid  2ir  radian  change.  Also,  as  will  be  shown  in  Section 
IV,  the  effect  of  loss  in  the  solid  medium  can  readily  be 
observed  in  and  estimated  from  the  reflectance  function 
R(6).  There  are  other  potential  areas  of  application  in¬ 
volving  film  thickness  measurement  and  the  thin-fllm 
acoustic  matching  between  two  media  of  vastly  different 
impedances,  which  will  also  be  explored  in  Section  IV. 

II.  Nonparaxial  Formulation  of  the  K(z)  Integral 

The  derivation  of  the  relation  between  V(z)  and  R(d)  in 
the  literature  [1],  [2],  [6],  [9]  can  be  classified  into  two 
general  categories:  those  involving  Fourier  spectral  de¬ 
composition  [1],  [2],  [6]  and  those  based  on  ray  optic 
models  [9J.  The  angular  spectrum  approach  is  more  ame¬ 
nable  to  inversion  since  Viz)  can  be  expressed  as  an  inte¬ 
gral  involving  /?(0),  and  the  integral  can  easily  be  manip¬ 
ulated  to  establish  a  Fourier  transform  relation  between 
F(z)  and  R(6).  The  ray  optic  approach  essentially  regards 
Viz)  as  the  interference  between  a  family  of  specularly 
reflected  rays  and  a  family  of  laterally  displaced  leaky  sur¬ 
face  acoustic  wave  rays.  It  is  not  complete  in  that  it  ne¬ 
glects  the  contributions  from  regions  near  the  longitudinal 
and  shear  critical  angles.  Also,  since  it  does  not  take  dif¬ 
fraction  into  account,  it  is  not  valid  near  the  focal  plane. 
In  addition,  the  ray  optic  formalism  in  the  form  usually 
given  is  incompatible  with  inversion  procedures  for  find¬ 
ing  R{6).  A  common  drawback  to  both  kinds  of  formula- 
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Fig.  1.  Schematic  of  spherical  focused  transducer  used  in  V(z)  measure¬ 
ments. 


tions  is  that  they  both  make  use  of  the  paraxial  approxi¬ 
mation. 

An  alternative  hybrid  formulation  of  the  relation  be¬ 
tween  K(z)  and  R(6)  will  be  given  in  the  following  treat¬ 
ment.  This  approach  combines  features  from  the  angular 
spectrum  as  well  as  the  ray  optic  models  and  is  valid  for 
all  z  including  the  focal  plane.  Moreover,  no  paraxial  as¬ 
sumption  has  to  be  made  in  its  development.  The  special 
case  of  a  focused  spherical  transducer  radiating  directly 
into  water  with  no  buffer  medium,  as  shown  in  Fig.  1,  is 
considered.  The  practical  significance  of  this  acoustic  lens 
configuration  will  be  discussed  in  Section  IV.  The  elimi¬ 
nation  of  the  buffer  rod  simplifies  the  derivation.  However, 
the  final  solution  is  valid  for  any  lens  configuration  that 
generates  a  converging  spherical  wave,  and  can  be  gen¬ 
eralized  for  use  with  other  systems. 

Using  the  notation  given  in  Fig.  1,  the  potential  or 
pressure  field  emitted  by  a  transducer  of  diameter  2  r0  is 
equivalent  to  that  produced  by  a  perfect  converging  spher¬ 
ical  wave  impinging  on  a  circular  aperture  of  radius  r0. 
All  acoustic  field  terms  will  be  assumed  to  vary  as  exp 
(jut).  Referring  to  the  aperture  plane  P  in  Fig.  1,  the  po¬ 
tential  field  distribution  in  cylindrical  coordinates  is  given 
by  the  relation 


<t>(r)  =  P{r) 


exp  [jkyJr2  +  f]] 


a) 


where  f\  =  V/o  -  d  or/,  =  f0  cos  0;  f0  is  the  focal 
length;  0O  is  the  maximum  half  angle  subtended  by  the 
aperture  at  the  focus;  and  k  is  the  propagation  constant  in 
the  liquid  medium.  P(r)  is  the  pupil  function  of  the  lens, 
which  in  the  most  general  case  takes  into  account  the  non- 
uniform  excitation  of  the  transducer  and  any  aberration 
effects.  In  the  simplest  case  for  an  ideal  lens 


P(r )  = 


0  <  r  <  r0 
r  <  r0 


(2) 


Making  use  of  the  circular  symmetry  of  the  field  distri¬ 
bution,  <f>(r)  can  be  decomposed  into  cylindrical  radial 
waves  of  wave  number  k,  by  applying  the  Hankel  trans¬ 
form 


A(kr)  =  ( 
Jo 


P(r) 


exp  [yWr  +  f\\ 


J0(krr)  r  dr.  (3) 


Here  Jfc2  =  ^  +  ifcj.  Since  the  transducer  diameter  is  nor¬ 
mally  many  wavelengths  wide  and  the  integrand  is 
weighted  by  r,  the  contribution  to  the  integral  by  Jo(krr) 
comes  mainly  from  regions  where  krr  is  large.  Therefore 
A(kr)  can  be  approximated  by  replacing  J0{k,r)  with  its 
large  argument  asymptotic  form 


exp 


Jofar)  - 


+  CXP  y  "  ?) 


yl2xkrr 


(4) 


Thus,  it  follows  that 


A(k.)  -  f  P(r)  ”PWL(ajL?iLO'!!lWI  rdr  (5) 
Jo 


s/2x  krr(i*  +  f]) 


where 


ti(r)  =  (kJr2  +f\  -  krr  +  ^ 
=  (k4? Tfi  +  k,r  -  £). 


h(r) 


The  method  of  stationary  phase  will  now  be  applied  to 
evaluate  (5)  [10].  Consider  the  phase  of  the  first  term  of 
the  integrand  i/',(r).  The  contribution  to  the  integral  comes 
from  the  region  near  where  ^J(r)  =  0 

*ito  *  -  *.  -  0 


or  where  r  =  r' 


kr 


/.  -  r/f 


(6) 


Note  that  r'  is  positive  in  this  case. 

Next,  consider  the  phase  of  the  second  term  of  the  in' 
tegrand  i/^fr).  Setting  V2(r)  to  zero  yields  the  result 


kr' 


~kr 


0) 


which  implies  that  r'  is  negative.  Since  in  a  cylindrical 
coordinate  system  r  is  always  positive,  there  is  no  contri¬ 
bution  to  A(kr)  from  the  second  term  of  the  asymptotic 
form  of  J0(kr  r) . 

Equation  (5)  therefore  becomes 


A(K)  = 


i  2*  exp  MMO]  kr , 

WHOI/  'J2xkrr'(r'1  +  /})  */' 


(8) 


where  r'  is  given  by  (6)  and 


Vti'(r’) 


kr'- 


(r':  +  f\)m 


k! 

A1/.' 


Substituting  (6)  into  (8)  gives  the  result 

A{kr)  =  P(kr)  k  exp  \jk.f\lkt 


(9) 
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where 


kz  =  Vfc2  - 


Mk- 0)  =  P(kr)  Ukt 


Mkn  z)  =  P(kr)  k  exp  [~jkzz)/kz. 


B(kr;  z)  =  R(kr)  A(kr\  z). 


V(Z)  -  S|| 


j  i<t>ru[  "  <« 

2 


Ur  = 


3z 

a*r 

dz ' 


=  - jk.A(kr ) 
W,)  =  jk.B(kr). 


Viz)  = 


/4(*,;  z)  Bikr\  z)  k.k,  dkr 


A(kr\  z)  A*(kr ;  z)  </*, 


Substituting  (12)  and  (13)  into  (17),  the  final  result  is 


(10) 


To  find  the  corresponding  angular  spectrum  at  a  plane 
of  distance  z  +  ft  from  the  aperture  plane  P,  each  spectral 
component  is  multiplied  by  the  propagation  factor  exp 
l ~jkt(z  +  /i)].  For  instance,  at  the  focal  plane  z  =  0 


Viz)  = 


! 


P\kr) 


(18) 


k,  dkr 


(ID 


which  simply  means  that  with  an  ideal  pupil  function  all 
angular  spectral  components  arrive  in  phase  at  the  focal 
point.  In  general,  for  any  arbitrary  z 


By  making  the  following  substitution  of  variables,  a 
comparison  can  be  made  between  (18)  and  Atalar’s  par* 
axial  angular  spectrum  formulation  of  the  V{z)  integral 
II].  Let 


(12) 


kjk  =  r/f0. 

Equation  (18)  can  be  written  in  the  form 


Suppose  the  plane  reflector  specimen  is  placed  at  z,  and 
the  reflectance  function  of  the  liquid-specimen  interface 
is  Rikr).  The  reflected  angular  spectrum  at  z  is 


P  V/o)  Rir/fo) 


Viz)  = 


exp  [-j2kzJ  1  -  jr/fp)1) 
Vl  ~  (r//0)* 


(19) 


r  dr 


P\r/f0) 


r  dr 


(13) 


Vl  -  ir/fo)2 


The  reciprocity  principle  of  Kino  and  Auld  given  in  the 
Appendix  is  employed  to  find  the  normalized  voltage  re¬ 
sponse  Viz)  of  the  focused  transducer  to  the  reflected  wave 
in  which  the  spectrum  is  given  by  (13).  Viz)  is  equal  to 
the  scattering  coefficient  5,, 


(14) 


where  uz  is  the  normal  displacement  at  the  liquid-solid 
interface,  and  the  integration  is  evaluated  over  the  same. 
The  superscripts  i  and  r  designate  the  incident  and  the 
reflected  fields  respectively,  and 

,  H‘ 


(20) 

which  differs  from  Atalar’s  solution  by  the  factor 
1/Vl  -  ir/fo)2  in  the  denominator  of  each  integrand.  In 
the  limit  of  the  paraxial  approximation  where  ir/f0)2  is  as¬ 
sumed  to  be  small,  this  additional  factor  tends  to  unity, 
and  (18)  becomes  the  expression  derived  by  Atalar. 

III.  Fourier  Transform  Relation  Between 
Viz)  AND  R(d) 

The  possibility  of  obtaining  R(d)  by  inverting  V(z)  has 
been  alluded  to  in  optical  microscopy  [11],  and  the  inver¬ 
sion  has  been  demonstrated  in  acoustic  microscopy  [12] 
based  on  Atalar’s  paraxial  Viz)  integral.  The  Fourier 
transform  relation  between  Viz)  and  Rid)  for  the  nonpar- 
axial  formulation  given  in  (18)  can  be  established  as  fol¬ 
lows.  Since  k%  =  Jr  -  k?r,  it  follows  that 


kzdk,  =  —  kr  dkr 
Equation  (18)  becomes 


(21) 


(15) 


Viz) 


-ci 


kcosdo 


P\kJRikz)  exp  ( -J2kzz)  dkz  (22) 


For  the  circularly  symmetric  incident  and  reflected  radial 
wave  components,  the  Hankel  transforms  of  the  uz  com¬ 
ponents  are 


where  C  is  the  normalization  factor  given  by 


-  *k cos0o 

-J* 


P2(kz)  dkz 


-i 


(23) 


(16) 


Noting  that  k,  —  k  cos  0,  where  k  =  2jt/X.  and  making 
the  following  change  of  variables 


By  applying  Parseval’s  Theorem  and  substituting  in  (16), 
(14)  can  be  written  in  the  form 


(17) 


z 

M  =  X 
t  =  2  cos  9 

we  obtain  the  result 

J  2  cos  0o 

°2 


(24) 


P2il)  Rit)  exp  i-j2irul)  dt.  (25) 
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In  short-hand  notation 


V{u)  =  ${P2(t)R(t)) 


(26) 


where  SF  {  }  denotes  the  Fourier  transform.  The  inversion 
to  obtain  R(t)  is  given  by 


rs-'im) 


R(t)  = 


0, 


P\t) 


2  >  t  >  2  cos  0O 
otherwise. 


(27) 


Note  that  R(t)  is  only  determined  for  the  angular  range  of 
insonification  covered  by  the  pupil  function  of  the  lens.  In 
general  F(r)  is  not  known  a  priori.  However,  /^(r)  can 
conveniently  be  calibrated  by  inverting  the  V(z)  function 
for  a  liquid-solid  interface  in  which  the  reflectance  has 
uniform  magnitude  and  phase  over  the  angular  extent  of 
the  lens  pupil  function.  In  particular,  if  a  perfectly  reflect¬ 
ing  surface  with  R(t)  =  1  is  used 


P\t)  =  SF"TOu)} 


(28) 


where  K0(u)  is  the  V(z)  of  the  “ideal”  reflector. 

This  inversion  approach  to  material  characterization  is 
different  in  that  it  requires  the  measurement  of  V(z)  in  both 
amplitude  and  phase  whereas  previous  work  [l]-[8]  uti¬ 
lized  only  the  amplitude  or  intensity  of  K(z).  It  is  generally 
not  trivial  to  conduct  high-accuracy  phase  measurement 
in  acoustic  microscopy.  A  special  phase-measurement  sys¬ 
tem  [13]  has  been  developed  to  enable  the  precise  acqui¬ 
sition  of  phase  data  and  the  subsequent  experimental  dem¬ 
onstration  of  the  inversion  technique. 

It  is  interesting  at  this  point  to  examine  the  range  res¬ 
olution  of  a  confocal  imaging  system  such  as  the  reflection 
acoustic  microscope  for  an  ideal  plane  reflector;  i.e.,  R(B) 
a  1.  Assuming  an  ideal  lens  pupil  function  of  the  form  in 
(2),  the  voltage  response  of  the  acoustic  microscope  as  a 
function  of  defocusing  distance  can  be  calculated  from 
(25)  to  be 

J  2  cos  0o 

exp  [~p.irut]  dt 


or 

K(«)  =  -  C  exp  [  —jxu(\  +  cos  0O)1 

•  sine  [2« (1  -  cos  0O)].  (29) 
The  3-dB  point  of  |K(m)|  is  given  by  the  relation 
2(1  -  cos  0O)  m,jb  =  0.457T 


or 


0.225A 

ZldB  ~ 

1  —  cos  0„ 

Therefore  the  3-dB  range  resolution  or  the  focal  depth  for 
a  plane  reflector  is 


^ZjdB 


0.45A 
1  -  cos  0o 


(30) 


IV.  Experimental  Results 

In  the  V(z)  measurement  a  special  acoustic  transducer 
is  used.  The  piezoelectric  element  is  a  spherical  shell  with 
an  air  backing  and  a  quarter-wave  matching  layer  on  the 
front  face.  The  transducer  radiates  directly  into  water,  and 
thus  forms  an  aberration-free  spherical  lens  that  focuses 
the  emitted  acoustic  beam  into  a  diffraction-limited  spot 
at  the  center  of  curvature  of  the  transducer.  This  acoustic 
lens  configuration  is  advantageous  for  V(z)  measurements 
because  the  elimination  of  the  buffer  rod  removes  the  sys¬ 
tematic  noise  due  to  reverberations  inside  the  buffer  rod 
that  tend  to  overlap  in  time  with  the  reflected  signal  from 
the  object.  With  such  a  construction  the  acoustic  lens  can 
be  translated  over  a  wide  range  of  z  to  collect  V(z)  data  in 
an  essentially  noise-free  environment.  The  maximum  limit 
to  the  vertical  movement  is  imposed  by  how  far  the  lens 
can  be  moved  towards  the  specimen.  The  problem  with 
multiple  echoes  in  the  water  path  between  the  transducer 
and  the  specimen  is  overcome  by  exciting  the  transducer 
with  a  short  tone  burst  and  time-gating  the  first  arrival 
return  signal. 

Two  transducers  of  this  type  with  FI  l  and  F/0.7  aper¬ 
tures,  respectively,  are  used  in  the  experiment.  The  radius 
of  curvature  or  focal  length  f0  is  16  mm  in  both  cases. 
From  geometric  considerations,  the  maximum  possible 
negative  defocusing  distance  is  given  by 

Zinin  *  “/I  =  -/oVl-l/^  (31) 

where  F  is  the  /-number  of  the  lens.  For  the  F/0.7  lens 
|Zminl  =  11.2  mm,  and  for  the  F/l  lens  |zmin|  =  13.9  mm. 
The  importance  of  zmin  will  be  discussed  in  Section  V. 

Both  transducers  have  a  fundamental  frequency  of  ap¬ 
proximately  3  MHz.  For  K(z)  measurements,  they  are  op¬ 
erated  at  the  third  harmonic  frequency  of  10.17  MHz  for 
the  F/0.7  lens  and  10.7  MHz  for  the  FI  1  lens.  Many  prac¬ 
tical  considerations  go  into  this  particular  choice  of  oper¬ 
ating  frequency  range.  The  inversion  formalism  is  most 
conveniently  applied  to  a  lossless  liquid  medium.  The  loss 
in  water  is  essentially  negligible  at  10  MHz  and  so  satis¬ 
fies  this  requirement.  Moreover,  the  wavelength  in  water 
of  about  150  ^m  is  sufficiently  large  to  render  any  irreg¬ 
ularity  and  instability  in  the  vertical  translation  mecha¬ 
nism  insignificant.  By  the  same  token  the  measurement  is 
less  subject  to  error  as  a  result  of  surface  finish  imperfec¬ 
tions  in  the  sample. 

A  lead  sample  is  used  to  calibrate  the  pupil  function  of 
the  lens.  With  a  longitudinal  velocity  =  2200  m/s,  the 
longitudinal  critical  angle  for  the  water-lead  interface  oc¬ 
curs  at  6  =  43°,  as  illustrated  in  the  theoretical  plot  in 
Fig.  2.  The  reflectance  function  is  essentially  flat  in  both 
amplitude  and  phase  from  0  =  0°  to  40°.  Fig.  3  shows  the 
magnitude  and  phase  of  the  K(c)  for  lead  from  ;  =  -  8 
mm  to  z  -  +8  mm.  The  phase  displayed  is  in  fact  that  of 
^(z)  with  the  linear  phase  term  2  kz  subtracted  out.  Fig.  4 
shows  P2(0)  for  the  F/0.7  lens  obtained  by  inverting  the 
experimental  K(z)  for  lead.  The  pupil  function  is  not  uni¬ 
form  in  either  magnitude  or  phase  especially  near  the  edge 
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Fig.  2.  Magnitude  and  phase  of  the  theoretical  reflectance  function  of  a 
water-lead  interface. 
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Fig.  3.  Experimental  f'(:)of  a  water/lead  interface  measured  with  a  10.17 
MHz  770.7  focused  transducer. 
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Fig.  4.  Calibration  of  the  pupil  function  of  the  10.17  MHz  F/0. 7  focused 
transducer. 

of  the  lens  where  an  accentuated  peak  occurs  around  6  = 
38°.  Premature  roll-off  of  the  lens  angular  response  is  also 
evident,  and  the  illumination  from  the  outer  periphery  of 
the  lens  is  weak.  The  experimental  P2(0)  for  the  Ft  1  lens 


Fig.  6.  Experimental  K(z)  of  water-fiised-silica  interface  at  10.17  MHz. 


shown  in  Fig.  5  exhibits  a  similar,  though  less  pro¬ 
nounced,  peak  at  6  =  23.5°.  These  anomalies  in  the  lens 
response  are  due  to  nonuniformity  in  the  thickness  of  the 
matching  layer  as  a  result  of  difficulties  encountered  in  the 
fabrication  process  [23]. 

It  should  be  pointed  out  that  this  simple  calibration  pro¬ 
cedure  can  be  utilized  as  an  alternative  method  for  char¬ 
acterizing  and  predicting  imaging  performance  of  focused 
systems.  Instead  of  measuring  the  point  spread  function 
at  the  focal  plane  as  is  the  case  in  most  conventional  tech¬ 
niques,  the  characterization  is  done  in  terms  of  the  pupil 
function  illumination.  This  approach  has  direct  relevance 
and  application  in  the  fabrication  of  spherical  acoustic 
lenses,  where  the  sphericity  of  the  lens  and  the  uniformity 
of  the  matching  layer  are  of  primary  concern. 

Fused  Silica 

The  experimental  V(z)  for  fused  silica  is  shown  in  Fig. 
6.  Fig.  7  displays  the  corresponding  /?(0),  which  is  ob¬ 
tained  by  inverting  the  F(z)  and  then  dividing  out  the  cal¬ 
ibrated  value  of  P2(9)  for  the  F/0. 7  lens.  The  theoretical 


Fig.  7.  Comparison  of  the  theoretical  and  experimental  reflectance  function  Fig.  9.  Comparison  of  the  theoretical  and  experimental  reflectance  function 
for  a  water-fused-silica  interface.  for  a  water-aluminum  interface. 

ALUMINUM  SILICON  NITRIDE 


Fig.  8.  Experimental  K(z)  of  a  water-aluminum  interface  at  10.17  MHz. 


Fig.  10.  Experimental  V(z)  of  a  water-silicon-nitride  interface  at  10.17 

MHz 


reflectance  is  also  shown  for  comparison.  A  value  of 
vw  =  I486  m/s  is  used  as  the  velocity  in  water  in  the 
inversion  algorithm.  There  is  good  general  agreement  be¬ 
tween  theory  and  experiment  in  both  magnitude  and  phase. 
The  experimental  shear  critical  angle  at  B  =  23.5°  and  the 
Rayleigh  critical  angle  at  6  =  25.85°  compare  extremely 
well  with  the  theoretical  values.  The  longitudinal  critical 
angle  is  not  reproduced  in  the  experimental  result  due  to 
the  problem  of  degradation  angular  resolution  caused  by 
the  spatial  truncation  of  the  V(z)  data.  This  difficulty  will 
be  elaborated  upon  in  Section  V. 

Aluminum 

Figs.  8  and  9  are  the  K(z)  and  R($)  for  aluminum.  Again 
the  agreement  between  theory  and  experiment  is  very 
good.  The  longitudinal  and  shear  critical  angles  show  up 
clearly  in  the  experimental  result,  in  this  case  at  $  =  13.64° 
and  6  -  28.64°,  respectively.  The  Rayleigh  critical  angle 
is  located  at  6  =  30.86°,  which  gives  a  phase  velocity  of 
2897  m/s  for  the  leaky  surface  acoustic  wave. 


Silicon  Nitride 

Figs.  10  and  11  show  the  F(z)  and  R(6)  for  hot-pressed 
silicon  nitride,  which  is  a  ceramic  material  with  high 
acoustic  velocity.  All  the  critical  angles  occur  at  low-in¬ 
cidence  angles,  and  only  the  Rayleigh  critical  angle  at 
0  =  15.0°  is  reproduced  in  the  experimental  phase  curve 
It  should  be  noted  that  the  dip  in  the  experimental  mag 
nitude  curve  for  silicon  nitride  and  smaller  dips  at  the  Ray 
leigh  critical  angle  in  Figs.  7  and  9  for  fused  quartz  and 
aluminum  are  only  an  artifact  introduced  by  the  spatial 
truncation  of  the  experimental  F(z)  data.  Similar  dips  can 
occur  as  a  result  of  shear  loss  in  the  solid,  an  effect  widely 
known  as  the  Rayleigh  critical-angle  phenomenon  [141, 
[151.  However,  the  acoustic  losses  in  all  three  materials  *. 
are  known  to  be  negligible  at  10  MHz,  therefore  the  ex 
perimcntal  dips  cannot  possibly  be  associated  with  loss 
As  will  be  shown  in  the  following  section  on  error  analy*  ^ 
sis,  this  experimental  artifact  is  well  understood  and  it*  ,^C 
some  cases  can  be  alleviated.  Hence,  the  inversion  tech- 
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Fig.  11.  Comparison  of  the  theoretical  and  experimental  reflectance  func¬ 
tion  for  the  water-silicon-nitride  interface. 


Fig.  13.  Comparison  of  the  theoretical  and  experimental  reflectance  rune 
tion  for  a  water-Plexiglas  interface. 


Fig.  12.  Experimental  P(z)  of  a  water-Plexiglas  interface  at  10.17  MHz. 

nique  does  not  preclude  the  characterization  of  such 
acoustic  properties  as  shear  loss,  temperature  dependence 
of  the  Rayleigh-wave  velocity,  and  various  contributing 
factors  to  the  Rayleigh  critical-angle  phenomenon. 

Plexiglas 

Figs.  12  and  13  are  the  V(z)  and  /?(0),  respectively,  for 
plexiglas  or  lucite,  a  low-acoustic  velocity  and  high  loss 
material.  This  example  illustrates  how  material  loss  can 
be  determined.  The  longitudinal  critical  angle  is  located 
at  0  =  32.7°,  which  corresponds  to  a  longitudinal  wave 
velocity  of  2750  m/s.  The  magnitude  of  the  reflectance 
peaks  at  0.7  rather  than  one,  as  it  would  be  for  a  lossless 
substrate.  The  amount  of  diminution  of  the  peak  level  de¬ 
pends  on  the  loss  for  the  longitudinal  mode.  The  theoret¬ 
ical  curve  is  fitted  to  the  experimental  one  by  varying  the 
longitudinal  loss  factor  QL.  The  QL  is  found  to  be  about 
50,  which  translates  into  an  attenuation  coefficient  of  232 
Np/m  or  20  dB/cm  at  10  MHz,  which  agrees  well  with 
other  published  values  [16]. 


Fig.  14.  Experimental  K(z)  of  a  water-teflon  interface  at  10.17  MHz. 
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Fig.  13.  Comparison  of  the  theoretical  and  experimental  reflectance  func¬ 
tion  for  a  water-teflon  interface. 


Teflon 

V(z)  and  R(0)  for  teflon  are  displayed  in  Figs.  14  and 
15.  As  expected  no  critical  angle  is  observed,  since  teflon 
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5ftm  GOLD  FILM  ON 
FUSED  SILICA 


Fig.  16.  Experimental  t'(z)  of  a  water  and  5-pm  gold  film  on  a  fused-silica 
interface  at  10.7  MHz. 
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Fig.  17.  Comparison  of  the  theoretical  and  experimental  reflectance  func¬ 
tion  for  a  water  and  5-pm  gold  film  on  a  fused-silica  interface. 


has  a  longitudinal  velocity  of  1400  m/s,  which  is  lower 
than  that  of  water. 

Thin-Film  Measurement 

The  following  example  demonstrates  the  potential  ap¬ 
plication  of  acoustic  microscopy  in  thin-film  character¬ 
ization  and  also  confirms  a  thin-film  matching  phenome¬ 
non  predicted  by  numerical  computation.  The  thin-film 
structure  used  in  this  example  is  a  5-jim-thick  gold  film 
deposited  on  a  fused-silica  substrate.  At  10  MHz  the  thick¬ 
ness  of  the  gold  film  corresponds  to  about  1.5  percent  of 
the  longitudinal  wavelength.  Figs.  16  and  17  are  the  plots 
of  K(z)  and  R(9),  respectively.  The  agreement  with  the 
theoretically  generated  R(6)  is  good.  R(6)  shows  a  null  at 
9  =  17.23°,  an  angle  between  the  longitudinal  and  shear 
critical  angles  of  the  fused  silica  substrate.  Physically  this 
means  that  the  incident  longitudinal  mode  in  water  cou¬ 
ples  very  strongly  into  a  bulk  propagating  mode  in  the 
fused  silica  substrate.  The  physics  of  this  coupling  is  not 
well  understood,  but  it  is  believed  that  the  longitudinal 


wave  in  water  converts  in  the  gold  film  to  a  leaky  Sezawa 
wave,  which  leaks  into  the  fused  silica  substrate  in  the 
form  of  a  propagating  shear  wave.  The  coupling  efficiency 
and  the  angle  of  incidence,  at  which  maximum  transmis¬ 
sion  occurs,  have  been  shown  to  be  functions  of  the  him 
thickness  by  numerical  analysis.  Thus  by  measuring  the 
reflectance,  one  can  get  information  about  the  film  thick¬ 
ness.  The  experimental  demonstration  of  this  thin-film 
phenomenon  also  opens  up  the  possibility  of  using  the  thin- 
film  matching  technique  in  the  design  of  acoustic  trans¬ 
ducers  to  improve  the  transmission  efficiency  between  two 
media  with  vastly  different  acoustic  impedances. 


V.  Error  Analysis 

The  introduction  of  experimental  artifacts  in  the  exper¬ 
imental  R(d)  due  to  the  spatial  truncation  of  the  V(z)  curve 
is  treated  in  this  section.  This  problem  arises  because  of 
the  finite  distance  over  which  V(z)  data  can  be  collected. 
Equation  (26)  shows  that  V(u)  and  [P2(f)  R(t)\  form  a  Fou- 
rier-transform  pair.  It  can  be  shown  by  Fourier-transform 
theory  that  since  [P2(r) /?(/)]  is  finite  in  the  t  domain,  V(u) 
has  to  be  infinite  in  the  u  domain.  Experimentally  V(z) 
can  only  be  obtained  for  some  finite-width  data  window. 
Therefore  the  actual  waveform  V'(z)  used  in  the  inversion 
is  a  truncated  version  of  K(«).  Thus 

V’(u)  =  V(u)  reel  (u/D)  (32) 

where  u  =  z/X,  and  D  is  the  width  of  the  data  window 
defined  in  terms  of  the  number  of  wavelengths  in  water. 
Noting  that  multiplication  in  the  u  domain  corresponds  to 
convolution  in  the  transform  domain  t,  the  reflectance 
function  obtained  by  inverting  V'(u)  is 

[P2(t)R(t)Y  =  [P2(t)  /?(!)]*  D  sine  ( Dt )  (33) 


where  *  denotes  convolution. 

The  effects  of  the  convolution  between  the  reflectance 
function  and  the  waveform  sine  (Dt)  are  twofold.  First,  the 
angular  resolution  in  the  t  domain  is  degraded.  Second, 
because  of  the  oscillatory  nature  of  sine  (Dt),  sharp  fea¬ 
tures  in  P2(t)R(t)  tend  to  generate  ripples  in  the  resulting 
[/^(f) /?(/)]'.  The  resolution  degradation  effect  can  be  es¬ 
timated  as  follows.  The  full  width  between  zeroes  of  the 
main  lobe  of  sine  (Dt)  is 

At  =  2  ID.  (34) 


The  At  is  essentially  the  transition  width  of  the  response 
to  a  sharp  step  in  P2(t)R(t)  and  hence  can  be  regarded  as 
the  worst  case  resolution  in  t.  Using  the  relation  t  =  2  cos 
0  given  in  (24),  we  get  the  following  expression  for  the 
angular  resolution  in  terms  of  the  angle  of  incidence  9 


A  e 


At 

2  sin  O' 


(35) 


Combining  (34)  and  (35),  the  resolution  in  0  in  terms  of 
the  width  of  the  data  window  D  is  found  to  be 


A9 


_ 1_ 

D  sin  9' 


(36) 
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Fig.  18.  Variation  of  angular  resolution  degradation  as  a  function  of  the 
incident  angle  and  the  normalized  data  window  width  D. 


SILICON  NITRDE 


Fig.  19.  Simulated  reflectance  function  of  a  water-silicon  nitride  interface 
obtained  by  inverting  spatially  truncation  F(z)  data. 


Equation  (36)  is  plotted  in  Fig.  18  with  D  as  a  param¬ 
eter.  The  data  window  D  used  in  the  experiment  is  about 
100.  This  implies  angular  resolution  of  about  3.3°  at  0  = 
10°,  1.7°  at  0  =  20°,  and  1.1°  at  0  =  30°.  The  experi¬ 
mental  reflectance  of  silicon  nitride  in  Fig.  11  is  a  patho¬ 
logical  example  of  the  degradation  effect  of  this  angular 
resolution.  The  Rayleigh  critical  angle  occurs  at  a  low 
value  of  0  =  15.0°  and  in  the  vicinity  of  this  angle,  the 
magnitude  of  the  reflectance  is  unity  but  the  phase  goes 
through  a  rapid  2k  radian  change  over  a  2°  angular  range. 
The  net  result  of  the  convolution  with  sine  (£>f)  is  a  sharp 
dip  in  the  magnitude  of  the  experimental  reflectance  and 
a  smoothing  of  the  phase  curve  near  the  Rayleigh  critical 
angle.  The  effect  of  the  convolution  is  much  less  serious 
at  angles  of  high  incidence  as  shown  in  Fig.  18  and  as 
evidenced  by  the  experimental  R(0)' s  of  fused  silica  and 
aluminum,  where  only  a  slight  dip  occurs  in  the  reflec¬ 
tance  magnitude. 

A  computer  simulation  has  been  carried  out  where  the 
theoretical  R(0 )  for  silicon  nitride  is  used  to  generate  the 
F(z),  which  is  truncated  and  then  inverted  in  an  identical 
manner  to  the  experimental  data.  The  result  is  shown  in 
Fig.  19.  The  simulated  R(0)  exhibits  exactly  the  same  be¬ 
havior  in  both  amplitude  and  phase  as  the  experimental 
one  in  Fig.  11. 

The  locations  of  the  critical  angles  are  generally  used  to 
determine  the  phase  velocities  of  propagating  modes, 
which  are  given  by 

v.  =  u^/sin  d.c 

where  vw  is  the  velocity  in  water;  the  asterisk  stands  for 
L  (longitudinal).  5  (shear),  or  R  (Rayleigh);  and  the  sub¬ 
script  c  denotes  critical  angle.  The  percentage  error  in  the 
estimation  of  v.  as  a  consequence  of  the  truncation  of 
F(z)  can  be  shown  to  be 


Av. 
v . 


1  cos  0.  _  I  cos  0 

2  sin  0.c  ^  *'  2D  sin2  0 


(37) 


Equation  (37)  is  plotted  in  Fig.  20  with  D  as  a  parameter. 


Fig.  20.  Phase  velocity  estimation  error  as  a  function  of  angle  of  incidence 
due  to  degradation  of  the  angular  resolution. 


The  error  decreases  drastically  with  increasing  incidence 
and  increasing  D.  For  D  =  100,  |  Au./o.|  is  four  percent 
at  0  =  20°  and  1.75  percent  at  0  =  30°.  One  should  bear 
in  mind  that  (37)  represents  the  worst-case  estimation,  and 
the  actual  error  could  well  be  substantially  smaller,  espe¬ 
cially  for  the  Rayleigh  critical  angle  0Rc.  The  determina¬ 
tion  of  0Rc  involves  locating  the  point,  where  the  phase  of 
R(0)  is  ir  radians.  Around  0Rc  the  magnitude  of  the  reflec¬ 
tance  function  for  a  lossless  material  is  constant  while  the 
phase  ¥  can  be  shown  to  have  the  form  (9) 


tan  ¥  =  2 


*r 


k  sin  0Rc 
a 


(38) 


where  a  is  the  leak  rate  of  the  Rayleigh  wave.  The  phase 
function  is  the  antisymmetric  in  kr  about  V  =  k  around 
the  point  kr  =  k  sin  0Rc.  It  is  similarly  antisymmetric  in  k. 
or  t  in  the  neighborhood  of  k,  =  k  cos  0Rc.  Since  the  sine 
function  with  which  (P2(/)/?(f)J  is  convolved  is  symmet¬ 
ric,  the  phase  function  V  around  the  Rayleigh  critical  an¬ 
gle  remains  essentially  unchanged  by  the  convolution. 
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provided  that  the  width  of  the  main  lobe  of  the  sine  func¬ 
tion  Ai  is  smaller  than  the  transition  width  of  ^  through 
the  Rayleigh  critical-angle  region.  For  example,  in  the  case 
of  fused  silica  and  aluminum,  the  respective  leaky  sur¬ 
face-acoustic-wave  velocities  calculated  from  the  experi¬ 
mental  phase  curve  are  almost  identical  to  those  predicted 
by  theory. 

The  second  problem  associated  with  the  generation  of 
ripples  can  be  partly  remedied,  though  with  some  further 
loss  in  angular  resolution,  by  applying  a  smooth  apodi- 
zation  function  for  the  V(z)  data  before  inversion  takes 
place.  The  apodization  function  used  in  processing  the 
K(z)  data  is 


a(u)  = 


sin2  (vu/D) 
n  +  sin2  (x u/D) 


(38) 


where  n  is  a  free  parameter  for  adjusting  the  tapering 
characteristics  of  the  apodization.  The  n  =  0.1  is  used  to 
process  the  experimental  data  shown  here.  The  resulting 
apodization  essentially  leaves  the  data  in  the  center  of  the 
window  unchanged,  but  it  behaves  much  like  Hamming 
weighting  at  the  edges  of  the  data  window. 

The  obvious  solution  to  both  the  resolution  degradation 
and  ripple  problems  is  to  increase  the  data  window  width 
D  of  V'(u).  Equation  (31)  shows  that  this  can  be  realized 
by  increasing  the  radius  of  curvature  /0  or  reducing  the 
aperture  size  or  both.  Also  D  can  be  increased  by  increas¬ 
ing  the  frequency  of  operation.  For  silicon  nitride,  R(6) 
can  be  reproduced  much  more  faithfully  if  AO  is  reduced 
to  i  ",  which  can  be  achieved  by  using  an  F/1.5  lens  with 
a  focal  length  of  32  mm  operating  at  20  MHz. 

Another  important  source  of  error  is  in  the  estimation 
of  the  velocity  of  water.  Since 

u  =  zJ\  =  zf/v 


then 


zf 


du  = - j  dv  =  -udv/v. 

v 


(39) 


Suppose  the  wrong  velocity  is  used  in  the  inversion.  In 
(25)  this  is  equivalent  to  changing  u  to 

u'  =  u  +  Am  =  u(l  -  A v/u).  (40) 

The  resulting  inversion  is  given  by 

[/*2(/) /?(/)]'  -  j  V(u)  exp  \j2md  1  -  A vlv)t]  du 
or 


that  the  shift  distortion  as  a  function  of  incidence  angle  is 
of  the  form 


|A0|  =  cot  0  Aviv  (42) 

which  shows  that  the  most  serious  errors  occur  at  low  an¬ 
gles  of  incidence. 

The  problem  of  determining  the  true  transform 
P2(0)R(0)  from  a  finite  segment  V'{z)  of  V(z)  is  a  common 
one  in  Fourier  analysis  and  spectral  estimation.  Various 
techniques  exist  in  the  literature  for  extrapolating  V'(z)  so 
that  a  more  accurate  determination  of  P2(0)R(0)  can  be 
made.  Since  P2(0)R(0)  is  a  bandlimited  function,  the 
maximum  spatial  frequency  being  limited  by  the  angular 
extent  of  the  pupil  function,  V(z)  is  analytic  in  the  entire 
z  axis  [17].  In  principle,  an  iterative  algorithm  proposed 
by  Papoulis  [18]  can  be  used  to  improve  the  accuracy  of 
the  estimation  of  P2(O)R(0). 

Conclusion 

We  have  demonstrated  that  the  reflectance  function  of  a 
liquid-solid  interface  can  be  determined  by  inverting  the 
corresponding  complex  V(z)  data  from  an  acoustic  micro¬ 
scope.  This  inversion  technique  represents  a  more  com¬ 
plete  approach  to  material  characterization  than  previous 
K(z)-related  work.  The  phase  velocities  of  the  various 
propagating  modes  in  the  solid  medium  can  be  obtained 
directly  from  the  reflectance  function.  The  effect  of  ma¬ 
terial  loss  can  also  be  observed  and  quantified.  Moreover, 
this  measurement  technique  provides  a  means  of  gauging 
imaging  performance  of  focused  systems  by  directly  mea¬ 
suring  the  pupil  function  illumination.  In  addition,  there 
are  useful  practical  applications  in  the  area  of  thin-film 
characterization,  and  many  interesting  possibilities  exist 
for  more  complex  structures,  such  as  multi-layered  films. 
The  nonparaxial  formulation  of  the  V(z)  integral  is  im¬ 
portant  in  that  it  lays  a  sound  theoretical  foundation  for 
the  inversion  measurement  technique.  The  excellent 
agreement  between  the  theoretically  and  experimentally 
obtained  reflectance  functions  further  supports  the  validity 
of  the  nonparaxial  theory. 

Although  this  work  has  been  carried  out  and  discussed 
in  the  context  of  acoustic  microscopy,  the  validity  and  ap¬ 
plicability  of  many  of  the  underlying  concepts  extend  to 
optical  microscopy  as  well.  Provided  one  can  accurately 
measure  (he  optica!  phase,  which  is  not  trivial  but  cer¬ 
tainly  realizable  [22],  the  inversion  algorithm  described 
here  can  be  used  to  obtain  the  optical  reflectance  function 
from  the  corresponding  complex  optical  l'(c)  function. 


[P2(f)/?(r)|’  =  P2P[(1  -  Au/u)/].  (41) 

Therefore  using  an  erroneous  v  results  in  a  stretched  re¬ 
flectance  function  in  t.  Since  the  velocity  of  water  has  a 
large  temperature  coefficient,  4  m/s  per  °C,  this  error  can 
be  significant.  It  may  cause  enough  misalignment  of  the 
inversions  of  the  sample  material  of  interest  and  the  cali¬ 
bration  lead  sample  to  have  a  serious  effect  on  the  deter¬ 
mination  of  R(0).  From  (35)  and  (41)  it  can  easily  be  shown 
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Appendix 

Auld  and  Kino  [  19]-[2I],  by  using  the  reciprocity  theo¬ 
rem.  were  able  to  determine  the  normalized  reflected  sig¬ 
nal  or  reflection  coefficient  s,,  from  an  object.  The  theory 
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for  longitudinal  waves  in  a  liquid  can  be  stated  in  the  form 
jo)  (  ( u'p  -  up')  n  dS 


where  the  integral  is  taken  over  the  surface  of  the  object, 
u  is  the  displacement,  and  p  is  the  pressure  fields  associ¬ 
ated  with  an  incident  wave  of  temporal  frequency  oi.  Here 
exp  (jwt)  time  dependence  is  assumed.  The  superscript  i 
denotes  the  incident  or  transmitted  wave  when  the  object 
is  not  present,  and  the  unsuperscripted  terms  denote  the 
total  field  at  the  obstacle.  The  parameter  P  is  the  power 
exciting  the  transducer  for  a  given  incident  signal  u\  p‘. 

The  total  fields  can  be  written  in  the  form 

u  -  u‘  +  ur  (A2) 

P  =  P‘  +  Pr  (A3) 

where  the  superscript  r  denotes  the  waves  reflected  from 
the  object.  Substituting  (A2)  and  (A3)  into  (Al)  yields 

ju  (  ( u‘pr  —  urp‘)n  dS 


If  the  object  is  a  semi-infinite  plane  normal  to  the  z  direc¬ 
tion,  j,,  can  be  written  in  the  normalized  form 


V(z)  = 


J  (u[pr  -  urzp‘)dS 

2  (  u\  p‘  dS 


where  the  asterisk  denotes  the  complex  conjugate.  Finally 
it  is  convenient  to  write  the  pressure  in  terms  of  the  po¬ 
tential.  For  a  liquid  with  u  =  V</>,  it  can  be  shown  that 
p  =  Hence  it  follows  that 


^(z)  = 


j  -  U^) 

2  j  dS 
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Note  that  for  a  perfect  plane  reflector  located  at  the  focal 

plane  ;  =  0  of  a  lens,  <t>r  =  <t>'  =  <£'*,  and  uz  =  -u‘z; 

therefore  F(z)  =  1. 
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|  ABSTRACT 

* 

f 

We  have  developed  a  more  complete  theory  for  the  V(z)  char¬ 
acteristic  of  an  acoustic  microscope.  This  theory  is  nonparaxial 
and  treats  the  effect  of  finite  ka  ,  where  a  is  the  radius  of 
the  aperture.  It  explains  well  the  asymmetry  of  the  V(z)  curve 
for  a  perfect  reflector  which  has  been  observed  experimentally  and 
has  been  difficult  to  explain  theoretically.  The  results  obtained 
are  also  of  importance  to  the  scanning  optical  microscope. 
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INTRODUCTION 

It  has  been  known  for  some  time  that  the  output  signal  from  a  scanning 
acoustic  microscope  operating  In  a  reflection  mode  varies  rapidly  with  the 
distance  of  the  focus  z  from  the  plane  of  the  reflecting  object.  The  vari¬ 
ation  of  the  parameter  V(z)  ,  the  voltage  output  as  a  function  of  z  ,  Is  the 
principal  contrast  mechanism  of  the  scanning  acoustic  microscope  when  used  to 
observe  solid  materials  in  a  reflection  mode. 

The  theory  of  this  mechanism  has  been  given  by  several  authors ,1“^  and 
developed  recently  in  a  nonparaxial  form  by  Liang  and  Kino.4  Most  of  the 
attention  has  been  devoted  to  the  Interaction  of  the  specularly-reflected 
waves  from  a  solid  substrate  with  the  Rayleigh  waves  that  propagate  along  It. 
Recently,  it  has  become  important  to  consider  the  same  type  of  V(z)  measure¬ 
ment  for  scanning  optical  microscopes.  Surprisingly  enough,  there  seems  to 
have  been  relatively  little  attention  paid  to  this  phenomenon  in  the  optical 
literature.5  In  this  case,  Rayleigh  waves  do  not  occur,  but  it  would  be 
useful  to  evaluate  the  resolution  in  the  range  direction  (the  z  direction) 
of  an  optical  microscope,  and  to  understand  the  phenomenon  fully,  so  that  the 
device  can  be  used  for  quantitative  measurements  of  height. 

The  equivalent  type  of  interaction  in  an  acoustic  microscope  occurs,  for 
Instance,  when  a  lead  sample  is  examined  in  a  water  medium.  In  this  case,  the 
Rayleigh  wave  has  too  low  a  velocity  to  be  excited.  We  and  other  authors  have 
found  that  in  such  experiments  the  V(z)  plots  do  not  obey  the  simple  theory 
that  has  been  developed.  In  particular,  the  V(z)  plots  are  asymmetric  about 
the  point  z  =  0  (focal  point  at  the  surface  of  the  object),  while  most  of 
the  theories,  with  the  exception  of  those  due  to  Bertoni,*’5  suggest  that  the 
function  should  be  symmetric. 
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We  have  found  that  there  are  two  basic  causes  of  the  asymmetry:  phase 
variations  at  the  lens  surface,  and  the  simplification  of  the  theory  to  use 
either  the  paraxial  approximation  or  asymptotic  assumptions  along  with  the  use 
of  the  method  of  stationary  phase.  This  Is  equivalent  to  assuming  that 
kr^  ■*>  •  ,  where  rg  Is  the  radius  of  the  lens  aperture,  and  k  *  2w/x  Is  the 
propagation  constant  of  the  wave  of  wavelength  x  .  In  this  paper,  we  show 
how  to  modify  the  method  of  stationary  phase  to  make  a . nonparaxial  calculation 
for  V(z)  and  take  account  of  the  finite  Integration  range,  which  is  limited 
by  the  edge  of  the  aperture.  The  results  obtained  agree  well  with  experiment. 


We  first  evaluate  the  cylindrical  spatial  frequency  spectrum  of  the  waves 
at  the  focal  plane  A(k  0)  ,  where  kr  is  the  radial  propagation  constant 
and  fj  the  distance  of  the  focus  from  the  plane  of  the  lens  aperture,  as 
Illustrated  in  Fig.  1.  We  then  obtain  V(z)  in  the  same  manner  as  In  refer¬ 
ence  4.  Second,  we  modify  the  analytical  expression  for  A(k  0)  by  taking 

*  * 

account  of  the  integration  limits,  then  use  a  modified  form  of  the  method  of 
stationary  phase,  and  compute  V(z)  as  before. 


NUIOICAL  COMPUTATION  OF  A(kr,0) 

We  consider  a  spherical  focused  transducer  and  use  scalar  theory  appro¬ 
priate  to  acoustic  waves  in  a  liquid.  This  approach  can  be  shown  also  to  be 
an  accurate  one  for  calculating  V(z)  for  an  optical  beam.  We  use  the  nota¬ 
tion  shown  in  Fig.  1,  taking  the  origin  z  =  0  at  the  focus.  In  the  aperture 
plane  z  =  -f^  ,  the  potential  field  distribution  in  cylindrical  coordinates 
is  given  by  Liang  et  al4  as: 


♦(r)  =  P(r) 


exp  [jk/r  +  f. 


✓r  +  f 
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where  =  /fg  -  ,  or  f^  =  fg  cos  0g  ,  fQ  is  the  focal  length,  rg  is  the 

radius  of  the  transducer,  9Q  is  the  maximum  half-angle  subtended  by  the  aper¬ 
ture  at  the  focus;  k  is  the  wave  number  In  water;  and  P(r)  is  the  pupil 
function  of  the  lens.  For  the  ideal  lens,  we  assume  that 


P(r) 


0  <  r  <  r 


0  •• 


(2) 


P(r)  =0  r  >  rQ 

Since  the  field  distribution  is  circularly  symmetric,  the  cylindrical  spatial 
frequency  spectrum  of  <j>(r)  at  the  aperture  plane  P  ,  denoted  by  A(kp  - 
fj)  ,  can  be  expressed  by  its  Hankel  transform: 

A(k 


/exp  [jk/r2  +  f2 

P(r)  - ^3= - 00(krr)  rdr 


u 

■/ 


I? 77f 


exp[jk/r2  +  f2] 


(3) 


o  "2  ♦  i 


J0(kpr)  rdr 


In  the  focal  plane,  it  follows  that 


A(kp,0)  =  A(kf,  -  fj)  exp  (-jkfj) 


(4) 


2  2  2 

where  kr  is  the  radial  wave  number,  and  k  =  k^  +  kz  ,  or  k^  =  k  sin  e  , 
kz  =  k  cos  e  ,  and  0  is  the  angle  of  the  wave  propagation  direction  to  the 
z  axis. 

The  parameter  A(k  0)  can  be  obtained  by  numerical  integration.  For 
n , 

this  integration,  the  steps  in  kf,  &(kr)  must  be  chosen  small  enough  in 
order  to  keep  the  phase  variation  of  the  integrand  reasonably  small  over  each 
step.  Therefore,  when  krg  is  very  large  (e.g.,  krg  >  1000),  the  computation 
time  tends  to  be  inordinately  long. 


The  amplitude  and  the  phase  of  the  numerically-computed  function  in  the 


focal  plane  A(kp,  0)  are  shown  in  Figs.  2a  and  2b  for  an  aperture  radius 
rQ  *  10.3  mm  ,  a  focal  length  fQ  =  16  mm  ,  frequency  f  =  10.7  Wz  ,  and  a 
wavelength  x  -  0.14  mm  with  water  as  the  propagation  medium;  this  corresponds 
to  a  value  of  krQ  =  462  .  For  comparison,  the  dashed  lines  show  the  corres¬ 
ponding  results  obtained  by  using  the  asymptotic  form  of  the  Bessel  function 
and  the  method  of  stationary  phase  to  integrate  Eq.  (3).4  This  latter  result 
can  be  expressed  analytically  as  follows: 


jk 

A(kr,0)  =  r 


A(kp  ,0)  =  0 


i  -  *2y 


k  <  k  sin  9n 
r  u 


k  >  k  sin  e. 
r  u 


This  simple  expression  predicts  an  average  value  of  A(kr,0)  both  in 
amplitude  and  phase,  but  does  not  predict  the  rapid  variations  (the  Fresnel 
ripple)  of  A(kr,0)  with  kp  .  The  simplifying  approximation  yields  an 
expression  for  V(z)  ,  which  is  symmetric  around  z  =  0  because  of  the 
constant  phase  with  kp  ,  and  does  not  always  agree  with  experiments. 


A  MORE  ACCURATE  ANALYTICAL  EXPRESSION  FOR  A(kr,0) 

When  krQ  >  1000,  the  numerical  integration  of  Eq.  (3)  is  very  time  con¬ 
suming.  Therefore,  It  is  advantageous  to  work  out  a  more  accurate  analytical 
expression  for  A(kp,0)  than  that  given  by  the  direct  application  of  the 
method  of  stationary  phase.  Following  the  analysis  by  Liang  et  al  we  use 
the  asymptotic  form  of  Jg(krr)  and  obtain  the  approximate  expression: 


exp  [jij>  1  (r)] 


WV-V-  •*.  A  .* 


0  /2nkrr(r2  +  f\) 
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where 


4»1(r)  =  k/r2  +  f2  -  k  r  +  — 
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Following  the  method  of  stationary  phase,  the  main  contribution  to  the  inte- 

I 

gral  in  Eq.  (6)  comes  from  the  region  ^(r)  =  0  ,  or  where 


r'  =  —  f1  »  f.  tan  9 
k  1  1 


Therefore,  ^(r)  can  be  expanded  in  the  region  near  r'  as: 


►j(r)  =  4»!(r * )  +  D(r  -  r')2 


where 


9  3 

1  k  cos  e 


1  1  k  co 

D  =  _.  (r.)  = - 

2  1  2  kf 


It  follows  that 


A(kr>  -  fj)  =  B(r') 


r'  exp  [j  *j(r')] 
'27Tkrr'(r’2  +  f2) 


where 


B(r')  =  J " exp  [jD(r  -  r')2dr 


In  reference  [4],  following  the  standard  approach  to  the  method  of  sta¬ 
tionary  phase,  the  integration  limits  were  taken  from  -•  to  »  ,  but  this  is 


not  a  valid  assumption  for  kr  0  and  k^/k  sin  eQ  .  Here  we  retain  the 
limits  of  integration  of  r  as  0  and  Tq  .  Although  the  expansion  of  Eq. 
(8)  is  not  accurate  over  the  whole  region,  it  is  still  adequate  since  the  ra¬ 
pid  phase  variation  of  the  integrand  causes  the  main  contribution  to  the  in- 

2 

tegral  still  to  come  from  the  region  where  D(r  -  r')  is  relatively  small. 
Thus,  we  can  evaluate  Eq.  (11)  with  finite  limits  and  write 
0 

B(r')  =  f  — exp  (jx2)  dx  + 

J  /D 

-✓Dr'  0 


/  i. 

J  /D 


2 

exp  (jx  )dx 


CC(X)  +  jS(X)  +  C( Y)  +  jS( Y)] 


(12) 


where  C(z)  and  S(z)  are  the  Fresnel  integrals  defined  as 


z 
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Therefore,  it  follows  that 


1 

A(k  ,  -  f.)  =  -  CC(X)  +  jS(X)  +  C(Y)  +  jS(Y)]  exp  j[k  f  +  ir/4]  (17) 

r  1  /I  cos  e  z  1 

s 

In  the  focal  plane,  the  complex  amplitude  of  the  cylindrical  wave  spectrum  is 


A(k  ,0)  =  — 1 - [C(X)  +  jS(X)  +  C(Y)  +  jS(Y)]  exp(— \ 

✓  2  cos  0  \4  / 


(18) 


Figures  3a  and  3b  show  the  amplitude  and  phase  of  A(kr,0)  in  the  focal  plane 
predicted  by  Eq.  (18)  for  the  same  case  shown  in  Figs.  2a  and  2b.  In  the  re¬ 
gion  near  k^/k  =  sin  ,  where  it  dominates  V(z)  ,  the  result  is  similar  to 
that  obtained  from  the  numerical  integrations  shown  in  Figs.  2a  and  2b. 


THE  RESULTS  FOR  V(z) 

Using  our  numerical  computed  and  analytical  results  for  A(kp,0) 

Shown  in  Figs.  2  and  3,  we  have  computed  V(z)  the  same  way  as  that  given  in 
[4].  In  Fig.  4a,  we  show  the  theoretically-predicted  value  of  V(z)  for 
lead,  obtained  from  the  numerical  calculation  of  A(kp,0)  .  Figure  4b  uses 
the  more  accurate  analytical  expression  for  A(kr>0)  [Eq.  (18)],  and  Fig.  4c 
uses  the  simple  expression  given  In  Eq.  (5),  as  derived  by  Liang.4  Fig.  4d 
shows  the  corrsponding  experimental  result.  It  is  clear  that  the  predictions 
given  in  this  paper  (Figs.  4a  and  4b)  agree  with  experiments  carried  out  by 
Liang  much  better  than  the  simpler  theory  given  in  ref.  [4]  (Fig.  4c). 

Both  Figs.  4a  and  4b  show  very  qood  agreement  for  the  locations  of  the 
maxima  and  minima  with  the  experimental  data. 

We  believe  that  the  remaining  disagreement  between  the  theoretical  and 
experimental  data  of  the  V(z)  comes  from  the  nonuniformity  of  the  spherical 
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transducer  (or  lens)  excitation,  amplitude  and  phase  weighting  of  the  excita¬ 
tion  of  the  transducer  or  lens.  To  show  the  effect,  we  tried  Gaussian  taper¬ 
ing  of  both  the  amplitude  and  the  phase  of  the  transducer  output  signal. 
Figures  5a  and  5b  show  the  tapering  functions  used  and  the  corresponding 
V(z)  for  the  same  case  as  in  Fig.  4a.  Our  reasonable  guess  for  tapering 
yields  results  in  Fig.  5b  which  agree  even  better  with  experiments. 

CONCLUSION 

This  paper  gives  a  better  theoretical  evaluation  of  V(z)  in  the  reflec 
tlon  acoustic  microscope  by  dealing  with  the  evaluation  of  A(kp,0)  more 
carefully.  We  have  postulated  that  the  remaining  disagreement  between  the 
theoretical  and  experimental  data  is  mainly  due  to  the  nonuniformity  of  the 
excitation  of  the  lens  (or  spherical  transducer)  and  shown  that  slight  and 
reasonable  changes  in  this  excitation  give  still  better  agreement  with  experi 
ments . 
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FIGURE  CAPTIONS 

Fig.  1.  Schematic  of  spherical  focused  transducer. 

Fig.  2.  Cylindrical  spatial  frequency  spectrum  of  wave  excited  by  a 
spherical  focused  transducer  (rg  =  10.3  mm  »  fg  =  16  mm  , 

A  =  .14  mm)  at  focal  plane. 

(a)  amplitude;  (b)  phase.  Solid  line:  Predicted  by  numerical 
computation.  Dashed  line:  Predicted  by  the  simple  expression  in 
[4]. 

Fig.  3.  Cylindrical  spatial  frequency  spectrum  of  wave  excited  by  the  same 
focused  transducer  as  In  Fig.  2,  at  focal  plane  predicted  by  the 
more  accurate  analytical  expression  (Eq.  18).  (a)  amplitude; 

(b)  phase. 

Fig.  4.  The  comparison  of  the  theoretical  V(z)  for  lead  by  using  A(k,0) 
shown  in  Figs.  2  and  3  and  the  experimental  results,  (a)  A(k,0) 

corresponds  to  the  solid  line  in  Fig.  2.  (b)  A(k,0)  corresponds  to 

that  In  Fig.  3.  (c)  A(k,0)  corresponds  to  the  dashed  line  in 

Fig.  2. 

Fig.  5.  (a)  Gaussian  tapering  function  of  transducer  output  signal. 

(b)  The  corresponding  V(z)  for  the  same  case  as  in  Fig.  4a. 
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